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SECTION  I 
INTRODUCTION 

During  the  last  few  years  numerical  methods  for  predicting 
the  flow  field  around  an  airfoil  have  advanced  tremendously.  The 
computation  effort  is,  however,  quite  large;  there  may,  therefore, 
still  be  room  for  methods  which  yield  results  with  acceptable 
accuracy  but  with  less  computational  work.  This  claim  is  made  for 
the  integral  equation  method,  originally  proposed  by  Oswatitsch 
(Ref.  1)  and  developed  further  by  Spreitzer,  Zierep,  Norstrud,  Hancock, 
and  Nixon  (Ref.  2  through  9).  The  claim  to  shorter  computing  time 
is  made  on  the  basis  of  numerical  experiments.  The  following 
observations,  which  probably  have  already  been  in  the  minds  of  the 
originators  can  serve  as  an  explanation.  The  integral  equation 
formulation  contains  terms  which  are  found  directly  from  the  Prandtl 
Glauert  approximation  and  additional  terms  which  originate  from 
further  (nonlinear)  compressibility  effects.  Frequently  these 
additional  terms  are  fairly  small;  it  may  be  permissible  to 
evaluate  them  with  relaxed  accuracy  requirements. 

The  present  article  analyzes  the  method  from  a  theoretical 
point  of  view.  The  author's  original  goal  was  to  identify  the 
feature  of  the  method  which  is  the  basis  for  its  success.  The 
explanation  which  he  found  is  nothing  more  than  a  paraphrase  in 
mathematical  terms  of  the  reason  just  given;  it  has  little  practical 
interest.  Further  developments  of  the  underlying  ideas,  namely 
the  interpretation  of  the  integral  equation  method  as  a  method  of 
weighted  residuals,  may  however  lead  to  improvements  of  the 
computational  procedure  and  this  should  make  the  present  study 
worthwhile. 


SECTION  II 


GENERAL  OBSERVATIONS  AND  A  SURVEY  OF  THE  ANALYSIS 

The  discussion  is  based  on  two  fairly  simply  observations. 

The  first  is  the  basic  motivation  for  the  definition  of  a  Green's 
function:  the  Green's  function  pertaining  to  some  linear 

differential  operation  (with  appropriate  homogeneous  boundary 
conditions)  constitutes  the  kernel  of  an  integral  operator  which 
inverts  the  original  differential  operator.  If  one  uses,  instead 
of  the  Green's  function  a  fundamental  solution,  then  one  obtains 
only  a  partial  inversion  of  the  differential  operator,  but  one 
still  accomplishes  a  reduction  of  the  dimensionality  of  the  problem 
by  one. 

With  the  exact  Green's  function  one  knows,  of  course,  the 
solution  of  a  given  problem  except  for  some  quadratures.  Let 
us  assume  nevertheless  that  one  carries  out  a  method  of  weighted 
residuals  with  the  Green's  function  as  weight  function.  Accordingly, 
one  assumes  some  form  of  the  approximate  solution  which  satisfies 
the  boundary  conditions.  It  will  depend  upon  a  number  (usually 
very  large)  of  unknown  parameters.  One  then  expresses  the  residuals 
within  the  flow  field  in  terms  of  these  parameters  and  forms  the 
integral  over  the  residuals  weighted  with  the  Green’s  function  for 
which  the  singular  point  lies  at  a  certain  value  of  (x,y) .  Applying 
Green's  formula  to  the  integral  over  the  whole  field  with  the 
point  (x,y)  excluded  by  a  small  circle,  one  finds  that  the  solution 
at  the  chosen  point  (x,y)  is  expressed  solely  in  terms  of  the 
prescribed  boundary  values.  One  thus  finds  the  exact  solution 
for  the  point  (x,y)  in  spite  of  the  fact  that  the  residuals  in  the 
field  do  not  vanish.  The  Green's  function  with  a  singular  point 
(x,y)  therefore  constitutes  a  weight  function  which  makes  the 
solution  at  the  point  (x,y)  insensitive  to  residuals  within  the 
flow  field.  Incidentally,  in  a  nonlinear  problem  one  must 
consider  the  differential  operator  which  arises  by  linearizing 
the  original  partial  differential  equation  for  the  vicinity  of  a 
certain  approximation ,  the  so  called  Frechet  derivative.  The 
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Green's  function  pertaining  to  the  Frechet  derivative  makes  the 
result  for  a  certain  point  (x,y)  insensitive  against  first  order 
terms  in  the  residual. 

If  only  an  approximation  to  the  Green's  function  is  avail¬ 
able,  then  the  residuals  within  the  flow  field  will  give  a 
contribution  to  the  expression  for  the  solution  at  the  point  (x,y). 
This  contribution  is,  of  course,  small  if  the  approximation  to 
the  Green's  function  is  close  and  then  it  is  permissible  to  evaluate 
the  effect  of  the  residuals  with  reduced  accuracy.  Actually,  this 
interpretation  is  not  much  more  than  the  argument  put  forth 
above.  An  approximate  Green's  function  is  available  only  in 
special  cases  (mainly  flows  which  are  symmetric  with  respect  to  the 
x  axis).  Further  work,  namely  the  solution  of  an  integral  equation 
is  necessary,  if  only  an  approximate  fundamental  solution  is 
available,  but  the  main  effect  of  the  use  of  such  weight  functions 
is  the  same. 

The  present  report  studies  a  second  observation  in  considerable 
detail.  We  shall  interpret  the  integral  equation  approach  to 
transonic  flow  as  a  method  of  weighted  residuals  in  which  the 
weight  functions  are  given  by  approximate  fundamental  solutions 
(and  in  favorable  cases  by  approximate  Green's  functions).  Each 
choice  of  the  singular  point  (x,y)  in  the  Green's  function  gives 
one  weight  function.  In  methods  of  weighted  residuals  one  has 
considerable  freedom  in  the  choice  of  weight  functions.  One  can 
therefore  replace  some  of  the  weight  functions  which  are  used 
in  the  integral  equation  method,  by  weight  functions  of  a 
different  character.  This  makes  it  possible  to  combine  the 
integral  equation  method  with  finite  element  or  finite  difference 
approaches.  This  is  advisable  in  regions  where  the  available 
approximation  to  the  ideal  fundamental  solution  is  poor. 

One  thus  is  led  to  alternative  versions  of  Nixon’s  extended 
integral  equation  method  and,  in  a  further  extension ,  to  an 
improvement  in  the  far  field  conditions  for  the  lower  transonic 
range.  This  last  development  ought  to  prove  useful  even  in  con¬ 
junction  with  conventional  finite  differences  approaches,  for 


it  makes  possible  a  reduction  of  the  size  of  the  portion  of  the 
flow  field  in  which  the  computation  is  carried  out. 


Section  III  of  the  report  is  a  compilation  of  different 
formulations  of  the  problem  in  terms  of  differential  equations. 

Section  IV  is  an  overview  of  different  integral  equation  formula¬ 
tions.  This  seems  to  be  desirable  not  only  because  we  look  at 
the  method  in  a  different  light,  but  also  because  most  publications 
treat  the  theoretical  side  rather  briefly  since  they  are  mainly 
oriented  torwards  the  numerical  aspects.  One  question  is  particularly 
intriguing:  does  the  integral  equation  method  fully  express  the 

shock  conditions,  or  it  is  necessary  to  express  part  of  the  shock  i 

conditions  separately.  Here  physical  intuition  cannot  be  considered 

, 

as  a  reliable  guide  because  of  the  global  character  of  the  integral 
equation  formulation.  One  may,  of  course,  resort  to  numerical 
experimentation,  and  this  has  been  done  in  the  past.  But  even  for 
those  who  are  primarily  oriented  toward  the  computational  practice, 
the  present  analysis  may  provide  reassuring  background  information. 

Only  the  important  results  are  compiled  in  Section  IV.  Details 
are  found  in  a  number  of  appendices.  To  someone  well  versed  in 
the  mathematical  techniques  used  in  this  context,  a  more  condensed 
version  would  probably  be  sufficient,  but,  at  least  some  of  the 
details  are  desirable  if  one  wants  to  show  the  subtle  differences 
in  the  assumptions  and  in  the  arguments  between  different 
formulations . 

Section  V  considers  the  conditions  that  are  to  be  imposed 
at  the  transition  from  a  subsonic  to  a  supersonic  flow  at  the 
sonic  line.  In  principle  the  condition  to  be  satisfied  has  local 

character  and  the  global  character  of  the  integral  equation  I 

formulation  is  not  too  well-suited  for  this  purpose.  An 
ingenious  device  has  been  introduced  by  Spreiter,  which  allows 

one  to  overcome  this  difficulty;  its  mathematical  meaning  will  ' 

be  explored.  The  modifications  of  the  method  which  are  discussed 
in  the  subsequent  sections  will  probably  make  the  use  of  this 
device  unnecessary,  but  it  is  of  interest  to  see  why  it  works. 


Section  VI  describes  the  original  form  of  the  integral 
equation  method, the  extension  proposed  by  Nixon,  and  the  modifi¬ 
cations  suggested  by  the  present  interpretation.  Section  VII 
examines  the  relation  between  the  integral  equation  formulation 
and  the  far  field  conditions. 

The  present  report  does  not  claim  to  give  a  complete  survey 
of  the  integral  equation  method.  Finer  points,  the  extension 
to  three  dimensional  problems,  and  to  problems  in  which  small 
harmonic  oscillations  are  superimposed  to  a  steady  field  will 
not  be  treated.  To  carry  over  the  ideas  of  the  present  report 
to  problems  of  this  kind  is  not  difficult.  In  the  practical 
realization  of  the  modification  suggested  here,  experiences  gained 
in  the  past  with  the  integral  equation  method  ought  to  prove 
useful . 


SECTION  III 

PARTIAL  DIFFERENTIAL  EQUATIONS  AND  SHOCK  CONDITIONS 


The  investigation  starts  with  the  potential  equation 
simplified  for  transonic  flow: 


(1  -  M2  -  <^-)T - +  ^ - =0  (1) 

x  xx  y  y  '  ; 

where  x  and  y  are  Cartesian  coordinates,"^  the  velocity  potential, 

M  the  free  stream  Mach  number,  and  k  a  constant  which  arises  by  a 

—  2—2 

development  of  1  -  ( y—  /a  )  for  the  vicinity  of  the  free  stream 
Mach  number;  a  is  the  free  stream  sound  velocity.  For  a  free  stream 
Mach  number  one,  <  is  equal  to  (y+  1)  where  y  is  the  ratio  of  the 
specific  heats.  The  equation  is  brought  into  a  standard  form  by 
means  of  the  Prandtl  Glauert  coordinate  distortion 

y  =  (1  -  M2)1/2y  (2) 

x  =  x 

and  by  setting 

<p  =  k(1  ~  (3) 


One  obtains 


(1 


■{.)<{>  + 

f  X  Txx 


4>  =  0 

yy 


(4) 


or 


3x^x 


rx 


)  + 


<J>  =  0  , 

Vyy 


(5) 


Eq.  (4)  shows  that  sonic  speed  occurs  for  =  1. 

In  simplifying  the  shock  conditions  for  transonic  flow,  one 
retains  the  requirement  that  the  momentum  equation  for  the 
direction  tangential  to  the  shock  be  satisfied.  This  requires 
that  the  potential  be  continuous  as  one  traverses  the  shock. 
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The  energy  equation  is  always  satisfied  if  one  computes 
the  pressures  from  the  Bernoulli  equation.  In  addition,  the 
requirement  of  conservation  of  mass  is  retained.  It  finds  its 
expression  in  the  potential  equation.  If  one  postulates  that 
across  the  shock,  the  potential  equation  's  satisfied  in  the 
weak  sense,  then  conservation  of  mass  is  guaranteed.  This  holds, 
of  course,  for  Eq.  (5)  as  well  as  for  the  original  potential  equation , 
Eq.  (1).  A  fourth  shock  condition,  preservation  of  momentum  in 
the  direction  normal  to  the  shock,  must  be  waved  if  one  works 
with  a  flow  potential.  The  errors  so  introduced  are  of  third 
order  in  the  shock  intensity  (expressed,  for  instance,  by  the 
jump  of  the  Mach  number  across  the  shock) .  This  is  in  keeping 
with  the  approximation  made  in  Eq.  (1). 

A  method  of  weighted  residuals  amounts  to  the  use  of  the 
concept  of  weak  equality.  The  shock  conditions, therefore, appear 
automatically,  if  one  solves  the  potential  equation  by  a  method 
of  weighted  residuals.  A  derivation  of  the  shock  conditions  using 
the  concept  of  generalized  functions  (for  which  the  concept  of 
weak  equality  plays  a  fundamental  role)  is  shown  in  Appendix  1. 

We  quote  the  result  in  the  following. 

Let  the  shock  be  given  by 

x  =  xg(y) 

and  let  [HJ*  the  jump  of  some  quantity,  H,  across  the  shock  from 
the  upstream  to  the  downstream  side 

[HI  _  =  H(xs(y)  -  i,y)  -  H(xg(y)  +  l,Y)  (6) 

r  >  0,  c  -*•  0 . 

Let  0  be  the  angle  of  the  normal  to  the  shock  with  the  x-axis. 

One  has 


dx 
_ s 

dy 


=  -tgB 
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(7) 


Conservation  of 
by 


mass  and  continuity  of 


<t>  2 

r  *  _JL_  1  + 

2 


-  [<f>  ]  + 

*  T  y  — 


t  across  the  shock  are  expressed 


(8) 


and 

=  o.  (9) 

An  equivalent  formulation  is  obtained  by  introducing 

4>x  =  u  (10) 

if  =  v. 

y 

(The  actual  velocity  components  can  be  obtained  from  u  and  v  by 
the  transformations  Eqs.  (2)  and  (3).) 

Then  one  has 

u  =  v  (11) 

y  x  ' 


and  from  Eq.  (5) 


5S  {U  -  J  )  +  fy  V  =  °* 

The  shock  condition  (8)  then  assumes  the  form 

[u  -  (1/2)u2]+  -  [v] \ (dxg/ dy )  =  0. 


(12) 


(13) 


An  expression  equivalent  to  Eq.  (9),  but  in  terms  of  u  and  v  is 

obtained  by  differentiating  Eq.  (9)  along  the  shock 

+  +  (14) 

(u]_(dxg/dy)  +  [v]_  =  0. 

One  might  be  tempted  to  derive  Eq.  (14)  from  Eq.  (11)  by  means 
of  the  concept  of  weak  equality,  but  Eq.  (11)  is  not  the  expression 
for  the  conservation  of  some  physical  quantity,  although  it  has 
divergence  form.  Such  a  procedure  would  therefore  not  be  physically 
justified,  although  it  gives  the  correct  result. 


:.,V  . 
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A  formulation  of  the  problem  solely  in  terms  of  u,  which 
has  originally  been  used  in  the  integral  equation  method,  is 
obtained  by  differentiating  Eq.  (5)  with  respect  to  x.  One  obtains 


„2  2  „2 

— 2^  (u  -  j  2  u  =  0.  (15) 

3x  3y 


In  the  approximation  used  in  all  of  these  derivations,  the 
boundary  conditions  at  the  profile  are  transferred  to  the  x  axis. 
Along  the  upper  and  lower  side  of  the  profile  the  values  of  v  are 
then  given  by  the  profile  shape.  In  the  formulation  in  terms  of  u 
the  values  of  uy  =  v^  are  known  at  the  profile. 

The  shock  conditions  (13)  are  expressed  solely  in  terms  of 
u  by  using  Eq.  (14).  Using,  in  addition,  Eq.  (7)  one  obtains 


[u] 


+ 


+  cos  2  3 


0. 


(16) 


To  express  Eq.  (14)  solely  in  terms  of  u,  we  differentiate 


it  along  the  shock  with  respect  to  y. 
operator 

dx 


3y 


_ s  3 

dy  37 


This  means  we  apply  the 


One  obtains 


3  ,  ,+  dXs 
3y  [u]-  dy 


3  +  dx 

+  37  (  dy 


,  d2  x 

S)2+[u]+  - * 


dy' 


3  r  i +  "XS 
37  fv]_  -  ° 


+  17  m! 


Here  Eqs. (11)  and  (12)  are  used  to  eliminate  v.  One  obtains 

3  S  .  c  ,  S,  O  r„  U  ,-T  ,  r..,-r  "S  _  *  Q  .  (17) 


,  dx 

2  5y  dF  + 


3  r  ,+  ,dxs,2  Dr  u2,+  ^  r  ,+  d x 

-37  lul-  «“dy)  '  37[U~2  ]-  +  [U]-  - 


dy 
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Eq.  (16)  and  (17)  give  the  complete  formulation  of  the  shock 
conditions.  In  practice  one  will  probably  avoid  the  formulation 
solely  in  terms  of  u  because  of  the  complexity  of  the  shock 
condition,  (Eq.  17) . 

If  a  u-field  satisfying  Eq.  (15)  (but  not  necessarily  Eqs. 

(16)  and  (17))  is  known,  then  one  can  always  construct  a  v-field 
by  the  use  of  Eq .  (11)  and  Eq .  (12).  One  determines,  along  some 

initial  line  x  =  const  which  intersects  all  lines  y  =  const,  the 
velocity  component  v  from  Eq.  (11)  by  an  integration  with  respect 
to  y.  It  is  best  to  choose  this  line  at  x  =  -  *.  There,  all 
derivatives  of  u  are  zero,  and  therefore  v  =  0.  One  then  determines 
v  by  integrating  Eq.  (11)  along  lines  y  =  const.  This  construction 
cannot  reach  points  which  lie  downstream  of  the  shock.  For  these 
points  one  must  again  use  Eq.  (12),  along  a  line  x  =  +°°  to  obtain 
initial  conditions  for  v.  Again,  one  obtains  v  =  0.  (Actually 
some  discussion  of  the  behavior  of  u^  at  large  values  of  y  is 
needed  to  arrive  at  this  result.) 

The  differential  equation  in  terms  of  u,  Eq.  (15),  arises 
from  the  equation  for  conservation  of  mass  by  a  differentiation 
with  respect  to  x.  Conservation  of  mass  is  therefore  guaranteed 
only  if  conservation  of  mass  is  satisfied  along  some  initial  line. 
The  construction  of  the  v-field  described  here  imposes  this  require¬ 
ment  and  in  this  manner  legitimizes  the  use  of  Eq.  (15).  Whenever 
a  u-field  has  been  computed,  one  can  assume  that  the  pertinent 
v-field  is  available. 
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SECTION  IV 

INTEGRAL  EQUATION  FORMULATIONS 


In  the  following  discussions  a  weight  function  denoted  by 
iJj  will  be  used.  The  following  expression  has  been  chosen: 

tH£-x,n-y)  =  d/2)  log((£-x)2  +  (n-y)2).  (18) 


The  umbral  variables  operating  in  future  integral  expressions  are 
the  variables  and  n  which  appear  in  this  equation.  They 
correspond  to  x  and  y.  Each  choice  of  the  point  (x,y)  defines 
one  weight  function.  The  weight  functions  satisfy 


0 


(19) 


except  at  the  point  (',n)  =  (x,y) .  These  weight  functions  may 

2 

appear  unconventional  because  they  change  sign  at  (£)-  x)  + 

2 

(n-y)  =  1  and  because  they  tend  to  infinity  if  the  argument  is 

zero  or  infinite.  Actually,  any  weight  function  w(£;,n)  can  be 
written  as 


w(f  n) 


2tt 


If 


Aw(Clf  nL)  ,r 


(.rx  - 


nx-  nJd^dr^ 


(20) 


where  A  denotes  the  Laplace  operator  and  f  and  nx  are  umbral 
variables  of  integration.  This  is  seen  by  a  standard  application 
of  Green's  theorem  to  the  equation 


Aw(f,1,n1)  =  f  (C1,n1) . 


One  forms 

n-^  'MC1~  f./hj^-nJdCj^dn-j^  =//f(/1,  n^)  <lj (^i“  f-*n1~n)df>1dn1. 

Assuming  that  w  and  its  first  derivative  vanish  outside  of  a 
certain  region,  one  obtains  for  the  left  hand  side  2uw(f,,n). 
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This  leads  directly  to  Eq.  (20)  .  The  set  of  weight  functions  p 
is  therefore  equivalent  to  conventional  weight  functions  of 
finite  support  (such  weight  functions  are  different  from  zero 
only  in  a  finite  region) . 

We  now  substitute  some  function  (again  denoted  by  <J>)  into 
the  potential  Equation  (5) .  If  the  differential  equation  is  not 
satisfied,  then  one  obtains  a  residual.  We  evaluate  the  integral 
over  the  residual  weighted  with  the  function  <jj .  The  resulting 
expression  is  denoted  by  Q(x,y).  The  arguments  (x,y)  in  Q 
refer  to  the  same  arguments  in  t p.  The  independent  variables 
x  and  y  in  Eq.  (5)  are  replaced  by  f  and  n.  Accordingly,  we 
form 


Q(x,y)  -  //  (<j>_-  4>r2/2))  +  ~  (<l>n)  ]<M£>x,n-y)d£,dn. 

After  certain  transformations  which  are  shown  in  detail  in 
Appendix  II,  one  obtains: 

L 

Q  (x,y)  =  2ir<fr(x,y)  -  /  M  £-x,-y)  (<|>  (  £+0) -$  U,- 0  ))d£ 

o  ( 22 ) 

o°  2 

+  f  xp  (  t-x,-y)  (<j>(  £,+0)~<M  f,-0))df  +(l/2)  ff  -^(f-x,n-y)4>£  U,n)dEd,i 

Here,  i2  denotes  the  entire  plane  punctured  at  the  point  (On) 

=  (x,y) .  The  profile  extends  along  the  x  axis  from  x  =  0  to  x  =  L. 

It  appears  as  a  cut  in  the  f , n  plane.  In  the  presence  of  a 
circulation  (denoted  by  )’)  the  cut  is  extended  along  the  wake 
(again  represented  by  the  x  axis) .  The  arguments  +  0  and  -  0 
refer  to  the  limiting  values  assumed  by  n  as  one  approaches  the 
cut . 

Eq.  (22)  does  not  contain  explicit  contributions  due  to  the 
shock;  the  shock  conditions  are  included  by  the  use  of  the  concept  of 
weak  equality  which  is  inherent  in  a  weighted  residual  formulation. 

A  failure  to  satisfy  the  shock  conditions  registers  in  the  expression 
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Q  in  the  same  manner  as  a  failure  to  satisfy,  in  a  certain  region, 
the  partial  differential  Eq.  (5) .  If  one  chooses  to  exclude  the 
shock  from  the  region ft  by  a  cut,  then  the  first  expression  is 
modified  in  the  following  manner.  We  denote  by  ft'  the  region  ft 
cut  along  the  shock  and  by  Q'  the  pertinent  function  Q.  The 
region  ft ' has  the  two  borders  of  the  cut  as  part  of  its  boundary. 

Then 

Q' (x,y)  =  Q(x,y)  +  ( f>g  ( n)  -x , n-y )  { [ <j>f  ( r,s  ( n)  ,  n )  ]  *  (-vs  (n)  , r]  + 

-  [<J>n  (ts  (n) ,  nl!(dcs/dn)  )dn  (23) 

-/  [<J>  (f._  (n)  ,n]  +  (i|v  (C.  (o)  -x,  n-y)  ( f  (n)  -x,  n-y)  (d£  /dn)  )dn. 
o  —  t,  s  ns  s 

The  integrals  are  extended  along  the  shocks  in  the  positive  n~ 
direction.  The  second  integral  vanishes  if  >f>  is  continuous 
along  the  cut.  The  first  integral  vanishes  if  the  shock  con¬ 
ditions  Eq. (8)  are  satisfied.  In  principle,  continuity  of  the 
potential  and  conservation  of  mass  are  separate  postulates,  if 
the  shock  is  excluded  by  a  cut. 

In  a  symmetric  flow  pattern,  one  has 

<p(C,+0)  =  <Mt,-0)  and  <^(£,+0)  =  -4>  (£,-0) 

and  the  second  integral  in  Eq.  (20)  vanishes.  <J>  at  the 
profile  is  given  in  any  case  by  the  boundary  conditions.  Re¬ 
writing  Eq.  (22)  one  then  obtains 

L 

Q  (x,y)  =  2 it <j>  (x,  y)  -2  /  ij;  ( f.-x  ,-y )  (  f,  ,+0)d.r. 

°  (24) 

+  (1/2)  //Tpr(C-x,  n-y)  (f.  ,n)df,dn. 

ft  "• 

The  same  equation  is  obtained  for  a  symmetric  flow  if  one  con¬ 
siders  only  the  upper  half  plane  and  uses  instead  of  ij>  the  Green's 
function  (whose  n  derivative  vanishes  at  the  axis) 
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j(!J'  ( ',-x,  n-y)  -  n+y)). 

Eq.  (24)  is  particularly  simple  because  the  ur  mown  function 
appears  only  in  the  double  integral. 

In  the  usual  integral  equation  method  the  direct  use  of 
the  velocity  components  u  =  and  v  =  4>  is  preferred.  This 
is  convenient  because  u  and  v  are  of  direct  physical  interest. 

These  equations  are  obtained  by  differentiating  Eq .  (22)  with 

respect  to  x  and  y.  Details  are  shown  in  Appendix  III  Eqs.  (A18) 
and  (A19).  One  obtains: 

L 

Q  (x,y)  =  2TTu(x,y)  +  /  (v  ( f, ,  +0) -v  ( t, , -0)  )  vr  (  f-x  , -y )  d\ 

o 

L 

+  /  (u(f :,+0)-uM  ,-0)  )»{■  (£-x,-y)df.  (25) 

o  n 

-  ( r/2)  u2  (x,y) -1/2  ffu2  (  f, ,  n)  tj’  ( £>x ,  ry-y)  dt'dn 

L 

Q  (x,y)  =  2ttv  (x,y)  +  /  ( v  (  £  ,  +0 ) -v  ( f, , -0)  ip  n ( ,‘-x  , -y )  df. 

1  o  1 

L 

-  /  (uf,+0)  -u(t,-0)  )  ipr  (C-x  ,-y)  dr  (26) 

0  ^ 

-(1/2)  //  u2  (T, /r))^^n(r,-x/n-y)df,dn. 

Q  1 

Here,  as  before,  u  refers  to  the  C»n  plane  punctured  at  the  point 
(t,n)  •-  (x,y)  .  The  first  integral  on  the  right  in  each  of  these 
expressions  is  known  because  v  is  given  as  boundary  condition. 

The  second  integrals  do  not  extend  over  the  wake.  The  last  terms 
arise  by  differentiation  of  the  double  integral  with  respect  to 
x  or  y  If  one  omits  these  terms  then  one  obtains  the  integral 
equations  of  the  linearized  problem;  and,  for  the  symmetric  flow, 
even  its  solution.  The  terms  and  ij'_  which  occur  in  the 

double  integralare  singular  at  (^,n)=  (x,y)  .  One  cannot  tell  in 
advance  how  significant  the  double  integral  is  in  comparison  to  the 
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terms  (u/2)u  (x,y) .  Nixon  has  found  by  numerical  experimentation 
that  a  useful  starting  approximation  for  an  iterative  solution  is 
obtained  by  omitting  the  double  integral. 

In  a  method  of  weighted  residuals  based  on  Eq.  (21),  one 
postulates  that  the  expression  Q(x,y)  vanishes  for  all  weight 
functions  (this  means  for  all  values  of  x  and  y) .  In  a  method  of 
weighted  residuals  based  on  Eq.  (25),  one  postulates 
that  Q  vanish  for  all  values  of  x  and  y.  At  the  end  of  Appendix 
II  it  is  shown  that  Q  is  continuous  at  the  shock.  Equivalence 
between  the  formulations  (Eqs.(2))and  (25))  is  therefore  obtained 
if  Q  =  0  along  some  initial  line  which  intersects  all  lines  y  = 
constant,  for  instance  a  line  at  x  =  -  «.  In  this  argument  Eq. 

(26)  does  not  play  a  role.  In  practical  computations  Eq.  (26) 
may  well  be  useful  for  the  determination  of  the  v  field,  but  in 
theory  the  irrotationality  condition  vx  =  u y  serves  equally 
well . 

Eq.  (25)  can  be  derived  from  a  number  of  different  (although 
related)  postulates.  The  flexibility  in  viewing  the  problems  so 
obtained  may  sometimes  be  useful.  In  Appendix  IV  it  is  shown 
that  Eq.  (25)  constitutes  the  weighted  residual  expression  (not 
only  the  x  derivative  of  such  an  expression)  if  one  applies  the 
weight  function  -i|>„  to  the  potential  equation. 

Appendix  V  uses  the  formulation  of  the  problem  in  terms  of 
the  velocity  components  u  and  v.  It  is  shown  that  the  expression 
(25)  is  obtained  by  applying  the  weight  function  -< (ir  to  the 

c, 

equation  of  conservation  of  mass,  Eq.  (12). 

Appendix  VI  shows  that  Eq.  (25)  is  obtained  by  applying  the 
weight  function  )(<  to  the  formulation  of  the  problem  in  terms  of 
u,  Eq.  (15).  These  results  are  not  unexpected.  However,  the 
manipulations  bywhich  they  are  substantiated  are  not  entirely 
straightforward. 

Eq.  (25)  gives  the  integral  equation  formulation  solely  in 
terms  of  u  and  v.  There  is,  however,  a  conceptual  difficulty. 

In  deriving  this  equation  one  has  assumed  that  u  and  v  are  partial 
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derivatives  of  a  potential.  The  existence  of  a  potential  is 
guaranteed  if  one  constructs  the  v  field  pertaining  to  a  given  u 
field  in  the  manner  shown  in  Section  III,  for  this  construction 
is  based  on  the  irrotationality  condition  u^  =  v^.  The  continuity 
of  the  potential  at  the  shock  is  expressed  by  the  jump  condition 
Eq.  (14) .  The  other  condition  to  be  observed  at  a  shock  is  con¬ 
servation  of  mass,  expressed  by  Eq.  (13).  One  will  ask  whether 
the  integral  equation  formulation  derived  from  Eq.  (25)  encompasses 
these  two  conditions.  One  might  be  inclined  to  dismiss  this 
question  by  the  observation  that  a  serious  error  in  the  formulation 
of  the  problem  as  the  omission  of  a  shock  condition  would  have 
been  discovered  during  many  practical  applications  made  with  the 
method.  The  theoretical  answer  is  given  at  the  end  of  Appendix  V. 

The  two  shock  conditions  are  indeed  satisfied,  provided  that  the 
expression  Eq.  (25)  vanishes  for  all  values  of  x  and  y. 

Appendix  VII  gives  even  more  specific  information.  It  shows 
that  the  shock  condition,  Eq.  (13) ,  which  expresses  conservation 
of  mass  is  already  satisfied  if  vanishes  for  points  (x,y) 
immediately  upstream  and  downstream  of  the  shock.  For  points 
(x,y)  for  which  Eq.  (26)  is  satisfied  immediately  upstream  and 
downstream  of  the  shock,  one  finds  that  a  certain  linear  com¬ 
bination  of  the  shock  condition,  Eqs .  (13)  and  (14)  is  satisfied. 

These  results  are  not  restricted  to  normal  shocks. 

There  is  no  need  to  introduce  separately  shock  conditions 

if  one  uses  Eq.  (25) .  To  the  contrary,  if  one  would  express  the 

shock  condition  a  second  time  by  an  independent  postulate  one 

might  introduce  ill  conditioning  into  the  formulation  of  the 

problem,  because  in  the  discretization  which  is  always  needed  in 

practical  work,  one  may  inadvertently  express  the  same  conditions 

★ 

in  two  forms  which  are  not  completely  identical. 

★ 

The  procedure  proposed  by  Nixon  in  Ref.  6  might  lead  to  this 
difficulty.  The  condition  Eq.  (29 c)  in  this  reference  is  a  con¬ 
sequence  of  conditions  Eq. (29a)  and  (29b)  and  can  therefore  not 
be  used.  In  a  later  publication,  Ref.  (  8  )r  an  equation  corresponding 
to  Eg.  (29c)  does  not  appear. 


SECTION  V 

CONDITIONS  AT  THE  SONIC  LINE 


In  nearly  all  approaches  to  transonic  flow  (finite  differences, 
finite  elements,  integral  equations  methods)  one  uses,  implicitly 
or  explicitly,  the  concept  of  weak  equality.  It  is  possible  that 
expressions  which  solve  a  partial  differential  equation  in  the 
weak  sense,  fail  to  satisfy  it  pointwise.  In  the  transonic 
case,  solutions  containing  an  expansion  shock  in  addition  to  a 
compression  shock  thus  become  possible.  One  must  introduce  into  the 
numerical  approach  the  requirement  that  expansion  shocks  are  not 
admitted.  If  one  has  an  approximation  to  the  flow  field  with  a 
smooth  transition  from  subsonic  to  supersonic  flow  and  one  studies 
perturbations  to  such  a  solution,  then  one  finds  that  there  are 
particular  solutions  which  contain  at  the  transition  point,  a 
logarithmic  singularity.  Only  in  the  weak  sense  can  expressions 
with  such  a  singularity  be  regarded  as  solutions.  Physically  such 
perturbations  are  not  admissible.  One  must  impose  the  condition 
that  at  the  sonic  line  the  partial  differential  equation  is 
satisfied  in  the  strong  sense.  In  the  two-dimensional  case  this 
is  simply  expressed  by  the  requirement  that 


It  is  shown  in  Ref.l)  that  the  sonic  operator  used  in  the  Murman 
Cole  procedure  can  be  interpreted  in  this  manner. 

In  the  framework  of  the  integral  equation  formulation,  the 
condition  of  a  smooth  transition  through  the  sonic  line  is 
expressed  in  very  ingenious  but  a  rather  indirect  manner  which 
has  been  introduced  by  Spreiter  (Ref. 2).  In  all  integral  equation 
methods  one  postulates,  of  course,  that  the  weighted  residual 
expression  (25)  be  zero.  This  expression  is  written  in  the  form 


i 


or 


u  -  u2/2  =  R 


2 

u  -  2u  =  -2R.  (28) 

This  is  a  quadratic  equation  in  which  the  right  hand  side  contains 
the  unknown  function  u.  In  the  iteration  procedure  used  to  solve 
this  problem  the  right  hand  side  is  evaluated  from  the  preceding 
iteration  step.  Solving  this  quadratic  equation,  one  obtains 

u  =  1  ±  /1-2R  .  (29) 

We  mentioned  before  that  at  the  sonic  line  u  =  1.  Therefore,  one  j 

must  have 

R  =  1/2.  j 

Besides,  in  order  for  u  to  be  real  upstream  and  downstream  of  the 
sonic  line,  R  must  have  a  maximum  at  the  sonic  point.  One  therefore 
has 

3  R/ 3  x  =0. 

(In  the  subsonic  region  one  then  uses  the  lower,  in  the  supersonic 
region  the  upper  sign  in  Ect.  (29)). 

If  the  maximum  of  R  is  smaller  than  1/2,  then  U  will 
never  reach  the  value  1,  and  the  transition  from  a  subsonic  to 
a  supersonic  flow  can  occur  only  in  the  form  of  an  expansion 
shock.  This  formulation  expresses  indeed  the  requirement  that  u 
reaches  the  value  1  as  one  approaches  the  supersonic  region  from 
the  subsonic  side  upstream. 

In  an  iteration  process  there  is  always  some  arbitrariness 
regarding  the  terms  relegated  to  the  right  hand  side.  The  form  . 

of  Eq.  (25)  might  suggest  that  one  uses  (instead  of  Eq.(28)),  an 
equation  of  the  form 

I 

u  -  u~/4  =  R. 
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Then 


u  =  2  *  /4-4 R . 


Here  one  finds  that  R  must  have  a  maximum  at  1  in  order  for  u  to 
pass  smoothly  through  the  value  2  along  a  line  y  =  constant.  Of 
course,  such  a  postulate  is  without  interest,  but  one  sees  that 
only  for  a  special  choice  of  the  iteration  procedure  will  Eq.  (28) 
yield  conditions  for  the  sonic  line. 

Spreiter  and  Oswatitsch,  and  zierep  and  Nixon  have  devised 
workable  procedures  to  incorporate  this  condition  into  the 
integral  equation  method.  Yet,  it  has  a  somewhat  artificial 
character;  a  combination  of  the  integral  equation  method  with 
finite  difference  or  elements  methods  which  becomes  possible,  if 
one  regards  the  finite  element  method  as  one  of  weighted  residuals, 
can  lead  to  a  more  direct  expression  of  the  conditions  at  the 
sonic  line,  based  on  Eq.  (27). 

Nixon  in  Ref.  8  has  given  a  formulation  similar  to  Eq.  (28) 
for  the  vicinity  of  the  shock.  Dividing  Eq.  (16)  by  [u]  +  ,  one 
finds 


4" 

u  +  u 


COS  f  -  1 


where  we  write  temporarily  that  u  and  u  for  the  values  of  u, 
respectively,  upstream  and  downstream  of  the  shock.  For  the  normal 
shock  ( >3=  0)  this  amounts  to  an  approximation  to  the  Prandtl 
relation 


where  w  and  w  denote  the  velocity  upstream  and  downstream  of  the 
shock  and  a*  is  the  sonic  velocity.  It  follows  that  in  a 
formulation  of  Eq.  (25)  in  the  form 


u  -  (  (cos  R)  / 2)  u  =  R. 
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R  is  continuous  through  the  shock  and  the  values  of  u  obtained  in 
each  iteration  step  by  solving  this  quadratic  equation  are 
connected  by  the  proper  shock  relations. 
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SECTION  VI 

THE  INTEGRAL  EQUATION  METHODS  AND  SOME  EXTENSIONS 


So  far  we  have  discussed  expressions  for  weighted  residuals 
formed  with  the  weight  functions  U .  If  one  has  an  expression  for 
I  (or  for  u  and  v)  for  which  the  original  differential  equation 
and  the  shock  conditions  are  satisfied,  then  the  weighted  residuals 
Q  (or  Q^)  will  vanish  for  all  values  of  (x,y).  One  will  remember 
that  (x,y)  are  the  parameters  which  characterize  the  individual 
weight  functions.  In  practice  the  condition  =  0  (Eq.  (25)) 
will  be  imposed  only  at  selected  lines  or  selected  points  of  the 
flow  field.  One  can  hope  that  the  desired  technical  information 
(the  pressure  distribution  over  the  profile)  can  be  obtained  by 
the  use  of  a  relatively  small  number  of  values  (x,y) ,  because  the 
weight  functions  used  here  (approximate  fundamental  solutions) 
have  the  property  of  making  the  expressi^  s  fairly  insensitive 
against  errors  in  the  residual. 

In  the  original  form  of  the  integral  equation  method  the 
condition  Qx  =  0  was  imposed  only  for  points  of  the  profile,  or 
perhaps  points  of  the  profile  and  some  points  of  the  x  axis  up¬ 
stream  and  downstream.  The  expression  Eq.  (25)  contains, 

besides  local  terms  and  expressions  which  can  be  computed  from 
the  boundary  conditions,  a  double  integral  which  contains  values 
of  u  from  the  flow  field.  This  is  the  only  term  in  which  a 
failure  to  satisfy  the  differential  equation  or  the  shock  conditions 
appears.  The  values  of  u  which  occur  in  this  double  integral  are 
now  approximated  in  terms  of  the  values  of  u  at  the  x  axis.  Here 
one  uses,  besides  the  values  of  u,  the  values  of  u^  given  at  the 
profile  by  the  boundary  conditions,  and,  if  one  wants  to  have  a 
better  approximation,  the  fact  that  at  the  x  axis  (as  everywhere 
else)  the  differential  equation  for  u  must  be  satisfied.  The 
expressions  which  one  chooses  must  vanish  sufficiently  quickly 
as  one  goes  tc  infinity.  From  a  practical  point  of  view  one 
prefers  to  choose  formulae  for  u  for  which  one  of  the  integrations 
needed  in  the  double  integral  can  be  carried  out  analytically. 
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Ultimately  the  double  integral  is  reduced  to  a  single 
integral  which  contains  as  unknowns , only  the  values  of  u  at  the 
points  of  the  x  axis  under  consideration.  Substituting  these 
epxressions  into  Eq.  (25)  and  for  unsymmetric  problems  also  in 
Eq.  (26)  one  obtains  one  or  two  integral  equations,  which  are 
subsequently  solved  by  iteration.  Under  favorable  circumstances 
the  contribution  of  the  double  integral  is  small  and  one  obtains 
in  this  manner  results  for  the  pressure  distribution  over  the 
profile  which  are  of  sufficient  accuracy  for  technical  purposes. 

A  very  attractive  feature  of  this  formulation  lies  in  the 
fact  that  the  possibility  of  the  occurrence  of  shocks  (in  this 
approximation  normal  shocks)  is  put  into  evidence,  if  one  solves 
Eq.  (25)  formally  as  a  quadratic  equation. 

In  cases  with  an  embedded  supersonic  region  u  is  definitely  not 
small.  At  the  sonic  line  one  has  for  instance  u  =  1.  A  more 
accurate  determination  of  u  is  therefore  desirable.  One  can,  for 
instance  introduce  for  u(x,y)  within  the  field  expression  which 
satisfy  at  x  =  0  the  partial  differential  equation  for  u  and  also 
the  first  few  of  the  derivatives  of  the  differential  equation  with 
respect  to  y.  This  procedure  is  limited  by  the  complexity  of 
the  expressions  for  u  obtained  in  this  manner. 

To  achieve  a  more  accurate  determination  of  u,  Nixon  and 
Hancock  impose  the  condition  =  0  not  only  along  the  x  axis 
but  also  along  some  lines  y  =  constant  within  the  field,  and  thus 
obtain  a  considerable  improvement  of  the  method. 

According  to  the  interpretation  given  in  the  present  report, 
one  then  uses  weight  functions  p  pertaining  to  inner  points 
(x,y)  of  the  flow  field.  The  author  believes  that  other  weight 
functions  may  serve  equally  well,  if  not  better.  To  obtain  good 
results  for  the  pressure  distribution  at  the  profile  by  the 
extended  integral  equation  method  it  is,  of  course,  necessary  that 
the  residuals  at  points  close  to  the  profile  be  small.  The  weight 
functions  1  take  into  account  residuals  throughout  the  flow  field, 
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the  local  residual  at  a  point  close  to  the  profile  may  fail  to 
vanish  although  Q  =  0  at  the  point  (x,y)  under  consideration . 

For  points  in  the  vicinity  of  the  profile  where  u  is  not  small, 
the  use  of  a  localized  weight  function  therefore  appears  to  be 
preferable  to  the  conditions  obtained  from  the  integral  equation 
formulation.  Following  this  idea  one  uses  finite  difference  or 
finite  element  formulations  for  a  region  which  extends  somewhat 
beyond  the  supersonic  region. 

Such  a  procedure  has  one  additional  advantage.  The 
expressions  Q  contain  contributions  from  all  points  of  the  flow 
field  at  which  u  is  to  be  determined.  In  contrast,  a  finite 
difference  formulation  connects  the  value  of  u  at  a  certain  point 
(x,y)  only  with  values  at  neighboring  points.  For  numerical 
purposes  this  is  a  significant  simplification. 

One  thus  arrives  at  a  hybrid  procedure.  For  those  parts  of 
the  flow  field  where  u  is  small,  one  may  take  advantage  of  the  fact, 
that  the  use  of  weight  functions  yields  equations  in  which  a 
dominant  portion  is  directly  expressed  by  linear  theory,  and  the 
remaining  nonlocal  (and  nonlinear)  terms  are  small.  For  these 
expressions  it  may  suffice  if  one  describes  the  values  of  u  which 
appear  in  these  correction  terms  by  means  of  a  very  small  number 
of  points  (x,y)  within  the  flow  field.  In  the  portion  of  the  flow 
field  where  u  is  not  small,  one  uses  a  finite  difference  or  a  finite 
element  procedure.  In  this  region  one  may  even  replace  the  simpli¬ 
fied  potential  equation  (Eq.  1)  by  a  full  potential  equation. 

Accordingly,  we  divide  the  flow  field  into  an  inner  region 
in  which  we  use  a  finite  difference  and  a  finite  element  procedure, 
and  an  outer  region  in  which  one  uses  the  integral  equation 
formulation  (Figure  10) .  The  inner  region  extends  beyond  the 
embedded  supersonic  region,  but  still  it  may  be  fairly  small.  The 
boundary  of  the  region  in  which  the  integral  equation  formulation 
is  applied  is  no  longer  confined  to  the  x  axis  along  this  boundary, 
the  values  of  v  as  well  as  the  values  of  u  are  unknown.  The 
derivations  shown  in  the  Appendix  include  cases  of  this  kind.  One 
obtains  the  following  expression,  Eq.  (A.17)  : 


(30) 


r 


Q  (x,y)  =  2  n  (u  (x,y)  -  (1/4)  u2  (x,y)  +  f  v('vcd0  +  y/dn) 

Dfi  n 

2 

-  (u-  ( 1/2)  U  )  ;  dn  -  ub  dr' 

w  r'  n  1 

-  (1/2  ,.  ru2df  dn  =0 

Here  the  region  o  is  the  part  of  the  flow  field  in  which 
the  integral  equation  formulation  is  applied,  punctured  at  the 
point  (x,y) .  If  (x,y)  is  a  point  of  the  contour  of  the  region 
in  which  this  formulation  is  applied,  then  one  must  carry  out 
a  limiting  process  particularly  in  the  second  integral  over  a  3. 

The  result  is  found  in  Eq.(A48).  It  is  probably  practical  not  to 
consider  the  boundary  o  in  Figure  11  as  a  line  which  separates 
different  regions  of  the  flow  field  and  for  which  a  matching  of 
the  solutions  must  be  carried  out,  but  rather  as  a  line  which 
separates  regions  in  which  the  residuals  are  weighted  in  a  different 
fashion.  For  points  of  this  boundary  one  may  apply  either  the 
weighting  of  the  residuals  by  means  of  the  integral  equation 
method  (in  this  case  one  would  apply  £q.  (A48)  or  one  may  use  the 

weighting  applied  in  finite  difference  methods.  Practical 
experience  must  show  which  procedure  is  more  practical.  One  has 
the  choice  of  including  in  the  integral  equation  formulation  the 
residuals  that  are  due  to  the  inner  region  (they  are  presumably 
small  because  one  postulates  that  their  local  averages  be  zero); 
then  the  boundary  between  the  two  regions  does  not  appear  in  the 
formulae,  and  the  region  0  covers  the  entire  flow  field.  Alter¬ 
natively  one  can  define  Q  to  be  only  the  outer  region.  The  latter 
form  is  probably  preferable. 

Because  of  the  nonlinearity  of  the  problem  iterations  cannot 
be  avoided.  But  even  in  extensions  of  the  method  where  one  deals 
with  linear  equations  one  will  probably  use  iterative  procedures  in 
the  outer  region  because  the  number  of  unknowns  is  fairly  large  and 
the  matrices  which  appear  are  not  sparse.  It  can  be  expected 
that  such  iterations  converge  well.  In  the  inner  region  one  will 
also  use  iterations  (for  instance  the  Murman  Cole  iteration)  at 
least  in  the  initial  stages  of  the  procedure  .  After  an  initial 
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approximation  has  been  attained  one  might  think  of  a  Newton 
procedure,  in  which  the  individual  linear  equations  are  solved  by 
direct  elimination.  So  far  these  important  practical  aspects  have 
not  been  explored. 


SECTION  VII 

INTEGRAL  EQUATION  FORMULATION  AND  FAR  FIELD  CONDITIONS 


In  practice  the  field  for  which  a  numerical  procedure  is 
actually  carried  out  is  finite.  The  effect  of  the  field  outside 
of  its  boundaries  can  be  approximately  taken  into  account  by 
means  of  the  integral  equation  formulation.  For  this  purpose  it 
is,  of  course,  necessary  to  express  the  values  of  u  in  the  distant 
field  (that  is  outside  of  the  boundaries  of  the  computed  field) 
in  terms  of  the  value  of  and  4>  or  alternatively  of  u  and  v, 
at  this  boundary.  The  situation  is  analogous  to  that  encountered 
in  the  original  integral  equation  method,  except  that  in  the 
distant  field  the  values  of  u  are  smaller  and  therefore  the 
approximation  is  less  critical.  Of  course,  in  this  case  the 
values  of  4  are  not  known  in  advance.  If  one  chooses  to  use 
finite  difference  or  finite  element  methods  throughout  the 
computed  part  of  the  flow  field,  then  the  integral  equation 
formulation  serves  solely  to  formulate  far  field  conditions. 

In  deriving  the  usual  far  field  conditions  one  assumes 
that  in  the  distant  field  the  linearized  potential  equation  is 
satisfied.  In  the  formulation  arising  from  the  integral  equation 
formulation  also  nonlinear  terms  in  the  flow  differential 
equation  for  the  distant  field  are  taken  into  account.  This  is 
an  improvement  which  ought  to  make  it  possible  to  reduce  the  size 
of  the  portion  of  the  field  in  which  the  actual  computations  are 
carried  out. 

Source  terms  appear  also  in  Klunker's  formulation  of  the 
far  field  conditions  (Ref.  12) ,  but  these  sources  lie  within  the 
computed  parts  of  the  flow  field,  while  in  the  distant  field 
it  is  assumed  that  the  linearized  potential  equation  is  sufficiently 
accurate . 

A  first  formulation  of  the  far  field  condition  is  obtained 
from  Eq.  (A6)  .  Here  \t  is  the  distant  field  and  its  inner 

contour.  This  equation  is  evaluated  for  a  =  1.  We  quote: 
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aiup(x,  y)  +  /  'j'  ( f.-x ,  n-y )  [  4>r  ( f, ,  n)  dn  -  <l>  (r,n)df] 
™2 


“  /  4>(€ ,  n)  ( '!* r  (r.-x,n-y)dn  -  (C-x, n-y) dfj 

3!i2  •  11 

-  F  (it  -  0  (x-L]L  ,y)  ) 

-(1/2)  ff  iJj (C-x, n-y)  (4>^  ( f, , n) df.dn  =  0 

ft  J 


(31) 


Provided  that  one  has  an  approximation  for  the  double  integral 
one  obtains  a  global  relation  between  and  its  normal  derivative 
along  3ft  .  The  expression  contains  a  term  due  to  the  circulation. 

A  form  which  comes  closer  to  the  far  field  conditions 
derived  (for  more  general  problems)  in  Ref.  13,  is  obtained 
if  one  chooses  the  point  (x,y)  inside  of  the  boundary  3f>  •  One 
then  must  set  >.k  =  0  in  the  last  equation.  The  expression  so 
obtained  agrees  with  Eq.  (39)  of  Ref . 13  except  for  the  con¬ 
tributions  of  the  nonlinear  terms  in  the  distant  field  ft,  and 
of  the  circulation.  The  latter  term  is  missing  because  the  author 
considers  in  this  reference  the  potential  after  ;he  contribution 
of  the  circulation  has  been  subtracted. 

To  obtain,  in  the  notation  of  the  present  article,  an  expression 
where  $  denotes  the  potential  with  the  circulation  included,  one 
writes  Eq.  (39)  of  Ref.  13  in  the  form 


/  (<t>(f,n)  -  l'[(l/2) 

3ft2 

n-y)df,}- 
0  (r. ,  n)  /2 it ]  ]  dn 


-  (<i  (f,n)  /2  tt  >  ] }  { (~ipr .(<.-x,  n-y)  dn 

/■J)  (f,-x, n-y)  {-  jy  t<m,n)  -  r  [  (1/2) 

+  [  <M  £ ,  n )  -  r[(i/2)  -  r.(r,fn)/2iil3dU 

r  P 


=  o 


Hence , 
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:t,n)  [-'I'f  (r,-x,n-y) dn  +t|>n(c,-x,n-y)dfj 


-  /  /  (f.-x,  n-y)  [-0  (r,n)dn  +  ^  (f,n)dr] 


-r{  /  ((1/2)  -  (0(/,n)/2ir))h^f  (C-x,  n-y)dn  (c-x,n-y)d£] 

-  n 

-/  \p  (f-x,n-y)  [-  ^  (1/2)  -  (0  (C, n)/2ir) ) d<c 

3n2 

+  ^  (  (1/2)  -  (6  (f,,n)/2n)  )d£]  }  =  0 

Now  one  observes  that9(f,n)  as  well  as  i|j(£-x,n-y)  satisfies 
Laplace's  equation  in  the  region  ft  .  (One  will  remember  that 
the  singular  point  (x,y)  lies  outside  of  the  region  ft.)  The  two 
integrals  with  the  factor  I’  extended  over  the  whole  contour  3ft 
therefore  give  zero.  This  contour  consists  of  the  large  circle 
3ft^,  the  contour  separating  the  distant  field  from  the  computed 
part  of  the  flow  field  Dft2 ,  and  the  two  borders  of  a  cut  along 
the  x  axis  It  is  shown  in  Appendix  II  that  the  integral 

3ft^  vanishes.  One  therefore  has 


Most  terms  in  the  integrands  vanish,  since  3ft.,  extends  along  the 


x  axis;  one  obtains 


+  r  /  [(1/2)-  (0  (r ,  n) /2tt  )  ]  ^  u-x,n-y)df.  =  r  /v  u-x,  n-y)dt 
3ft3  "  rJ;L  1 


=  -r  ( 7T—  0  (x-L.  ,y)) 


Hence,  from  Eq.  (32) 
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f  <l>  ( i »  n  5  l -f  r  ( >'  -x )  ,  n-y )  dn  (C-x, n-y) di] 

3fi2 

-  /  v(C-x,n-y)(  -fr(Cfn)dn  +  <fn(£-x,n-y )dC] 

3fi2 

-I’  (n-p  (x  -  L^,y)  )  =0. 

This  is  Eq.  (A6)  for  u  =  0  but  without  the  source  terms. 

Source  terms  in  the  Laplace  equation  can  be  included  also  in 
the  procedure  of  Ref.  13.  One  may  do  this  by  inspection,  or 
by  the  following  rationale.  One  divides  the  distant  field  into 
a  very  distant  field  in  which  no  source  terms  are  presented  and  in 
a  moderately  distant  field  in  which  sources  play  a  role  (Fig. 11) 

The  inner  boundary  of  the  very  distant  field  consists  of  a  portion 
denoted  by  Ofig  and  a  cut  along  the  x  axis  extending  from  to 
infinity.  Eq.  (33)  holds  along  '.rather  than  S^)  which  is, 

of  course,  the  outer  boundary  of  the  moderately  distant  field. 

Now  one  applies  Green's  theorem  to  Eq.  (5)  with  the  weight  function 
!j>  ( t-x,  n~y)  ;  (X,y)  outside  the  very  distant  and  moderately  distant 

fields.  Using  Eq.  (33),  one  arrives  at  an  expression  in  which 
only  the  contribution  of  3^2  occurs.  This  is  Eq.  (31)  with  a=  0. 

If  one  solves  for  the  inner  field  iteratively,  by  means  of 
a  Murman  Cole  iteration,  say,  then  one  needs  local  rather  than 
global  boundary  conditions  along  To  obtain  such  conditions 

we  follow  the  argument  used  in  Ref.  13. 

Assume  that  we  possess  in  the  distant  field  a  solution  of 

ib  +  <f>  -ci(x,y)=o 

*xx  yy  !  ' 1  u 


where  as  always 


q  =  (1/2)  4x2). 

We  define  an  extension  of  this  distant  field  throughout  the 

inner  field  which  satisfies  if  +  <j»  =0.  (Because  of  the 

xx  yy 
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presence  of  the  circulation  one  needs  a  cut  along  the  x  axis 
also  in  the  inner  field.)  The  values  of  4>  match  at  the  boundary 
between  the  inner  and  the  outer  field.  There  will  be  a  jump  of 
the  normal  derivative  at  the  boundary  between  the  two  fields.  This 
jump  can  be  represented  by  a  source  distribution  with  the  unknown 
intensity  f(s)  (or  f(u))  where  s  (or  o)  measures  the  distance 
along  the  boundary  line.  One  thus  obtains 

r  ( x , y )  =  (2n)-1{  /  f(o)  v  (Mo)  -x,  n(o)-y)  da 

Zil  (34) 

+  .r/q(£,n)<M£  -x,  n-y)  df.dn  }-t  r  (  (1/2)  -  0  (x(s)  ,y  (s))/2tt) 

Here  ‘71  denotes  the  distant  field  and  30^  its  inner  boundary,  the 
cut  along  the  positive  x  axis  is  not  part  of  The  form  in 

which  the  circulation  enters  is  taken  from  Appendix  II.  It  does  not 
matter  whether  one  writes  0(x,y)  or  G(x-L,y)  as  long  as  the  point 
where  0  is  singular,  namely  x  =  L,y=0,  lies  outside  of  the  distant 
field  (inside  the  inner  field)/ for  in  the  distant  field,  the 
difference  0(x-L,y)  -  0(x,y)  can  be  represented  by 

f  f(a)  a  ( f,  ( a )  - x ,  n  (o)  -y)  do. 

A  change  of  L  therefore  amounts  to  a  change  of  the  function  f(s). 

The  last  expression  and  its  normal  derivative  are  evaluated 

for  points  x(s)  and  y(s)  on  3fi,  the  boundary  between  the  distant 

field  and  the  inner  field  for  which  the  differential  equation  is 

integrated  numerically.  The  function  f(s)  must  be  determined  in 

such  a  manner  that  the  value  of  <b  and  $  from  the  outer  field 

n 

and  from  the  computed  part  of  the  flow  field  match.  The  last 
equation  can  be  used  directly  to  determine  41.  The  evaluation 
of  .“h  is  shown  in  Appendix  VIII.  We  repeat  the  definitions  and  the 
result.  Let  f(s)  be  the  angle  of  the  contour  3V  of  the  distant 
field  with  the  x  axis.  Then  one  obtains,  according  to  Eg.  (A52) 
and  (A53) , 


f 
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4>n  (x  (s)  ,y  (s))  =  (  f  (s)  /2) 

+  (2tt)  1  {  /  -f  (a)  ip  (a  ,x  (s)  ,y  (s)  )  cos  (3  (o) -0  (s) )  do 
F)Q  n 

+  /  -f  (a)  ip  (a,x(s)  ,y(s)  )sin(B(o)  -  3(s))da 

dQ 


-cos6(s)  ff  q(4,  n)  >Pr(C-x(s)  ,  n-y  (s)dSdn 


+sin3(s)  ff  q(£,nk‘  (C-x(s)  ,  n-y(s)d£dn 

Q  n 


-(  T/r  (s)  )  sin  (6(s)  -0(x(s),y(s))} 


(35) 


<t>  =(2tt)  1{  -  f  f  (  a)  ip  (o,  x  (s)  ,y  (s)  sin  (  0(  o) -0  (s)  )  da 

^  'N  n  ^ 


+  /  f  (o)  i|  (o,x(s)  ,y  (s)  cos  (P(o) -B(s)  )  da 


-s  in  3  (s)  ff  q  (  r>,  n )  Pn  ( C~x,  n-y)d^dn 

n  0 

+cos  3  (s)  //  q(f,,n)ip_  ( :>x,  n-y)dr,dn 

n  n 

-  (F/r)cos ( 6(a)  — 0 (x (s) ,y (s) ) } 


where 

r(s)  =  (x(s)2  +  y(s)2)1//2 

0(x,y)  =  arctg  (y/x)  ,  0  £  0<2*  (3?) 

Vn  (cf,x,y)  =ipf  (f  (o)  -x,  n(o)  -y)  cos  6  (a)  +ip  (f.  (o)  -xf  n  (o-y)  sin0  (a) 

<1^  (o,x,y)  =ipc  (f,  (o)  -x,  n(o)  -y)  sin3(o)  -tp  (f,  (o)  -x,  n  (o)  -y)  cos0  (o) 

For  special  choices  of  the  contour  (3=  0,  or  6=  n/2)  this 

expression  will  simplify  greatly.  Some  remarks  regarding  the 

2 

approximation  of  q  (in  essence  of  u  )  which  occur  in  the  double 
integrals  of  Eqs.  (31)  ,  (34)  ,(35)  ,  and  (36)  will  be  mac’e  in  the 
next  section. 

These  formulae  will  be  used  in  the  following  manner. 

Assume  that  one  has  an  approximation  to  the  function  f  (s) .  One 


.  -  -  • 
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I 


then  computes  the  values  of  f  and  ij>  from  Eq .  (34)  and  (35) . 

In  one  omits  the  contribution  of  the  circulation.  Using  the 

value  of  t'  as  boundary  condition  at  the  outer  edge  of  the 
computed  flow  field  one  computes  the  inner  field,  for  instance 
by  the  Murman  Cole  iteration.  It  solves  Eq.  (5)  with  the 
boundary  condition  assigned  at  the  profile  and  determines  the 
circulation  i'  so  that  the  Kutta  condition  is  satisfied.  The 
expression  for  f  at  the  outer  boundary  must  be  modified  if  the 
circulation  is  modified.  From  the  computation  in  the  inner 
region  one  extracts  the  values  of  f  at  its  outer  boundary, 
but  without  the  contribution  of  the  circulation.  Let  J> 

n ,  inner 

and  :  .  be  the  value  of  f  (with  circulation  omitted)  for 

n, outer  n 

the  respective  fields.  A  correction  Af(s)  to  f(s),  which  can  be 
expected  to  lead  to  a  convergent  procedures,  is  then  given  by 


A  f  ( s ) 


n , outer 


(s)  - 


n ,  inner 


(s) 


(38) 


This  discussion  refers  to  computations  in  terms  of  the 
potential  f •  The  procedure  for  a  procedure  in  terms  of  the 
velocity  components  u  and  v  is  very  similar.  The  counterpart  to 
Eq.  (31)  is  Eq. (A48),  if  (x,y)  is  taken  on  the  boundary  between 
the  inner  and  outer  region,  and  Eq.  (A43)  with  a  =  0  and  r=0  if 
(x,y)  is  taken  in  the  inner  region.  We  quote: 


u:  „  ( -x.  n-y)  [u  ( -  ,  g)  -  (l/2)u2  g)  ]dg  -  v(; ,  n)dt"  ] 


+  '  1  (~-x,n_y)  [  u  ( r  ,  n)  d>-  +  v(,-.,ti)dr|)  (39) 

,>;?  ' n 

2 

-(1/2)  //  .  .  1  -x,  ,')u2  (  >  ,  n)d;dn  =  0 


(x,y)  in  the  computed  region,  and 


■A 
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:i[u(x(s),y(s) 


-  (1/2)  sin2B (s) u2 (x  (s)  ,y  (s) )  ] 

[  u(  (o)  ,  n  (  o)  )  -  (  1/2)  COS  2  |3  (S)  U2  (  (0  )  ,  1 1  (  ■))!,  r  (  f  ('-)  -  X  (s)  ,  r  (  •  ) 

2 

-  y(s))do  +  /  [  v  (f. ,  n)  -  (1/2)  sin  6  (s )  cos/  (s )  u2  (  i  (a ) ,  n  (s )  )  j  (40) 

'2 

Vx.(Uoix(o))do  -  (1/2)  //v  _  (;',-x,n-y)  u2  (/,  n)  dn  =  o 

^  (d ) 

’  in  the  punctured  outer  region,  (x,y)  on  the  countour 
Notice  that  in  these  equations  a  contribution  due  to  the  circula¬ 
tion  does  not  appear.  (This  holds  only  for  the  steady  problem; 
in  the  unsteady  problem  u  is  not  continuous  across  the  wake.) 

An  alternative  formulation,  in  which  the  distant  field  is 
represented  by  a  source  distribution  in  its  interior  and  at  its 
inner  boundary  is  possible,  also  if  one  works  with  u  and  v. 

One  computes  for  an  approximate  function  f(s)  and  an  approximate 
circulation  the  values  of  and  at  the  inner  boundary  of  the 
distant  field  from  Eqs .  (35)  and  (36)  and  then  expresses  <•■  in 

terms  of  u  and  v.  The  inner  field  to  which  the  finite  difference 
procedure  is  applied  is  computed  with  the  values  of  4>  (expressed 
in  u  and  v)  as  boundary  condition.  At  the  outer  boundary  again  an 
adjustment  of  the  circulation  is  needed.  The  subsequent  procedure 
is  the  same  as  before. 


33 


SECTION  VIII 

APPROXIMATION  OF  u  IN  THE  DISTANT  FIELD 

In  application  of  the  far  field  conditions  an  approximation 
2 

for  q ( x , y )  -  (1/2) • (u  )/  x  is  needed  for  the  distant  field.  One 
wants  the  computational  effort  for  deriving  such  an  approximation 
to  be  small.  On  the  other  hand  the  purpose  of  the  far  field 
conditions  is  to  reduce  the  size  of  the  portion  of  the  flow  field 
in  which  accurate  computations  must  be  carried  out,  and  for  this 
purpose  a  good  approximation  for  q  in  the  distant  field  is  needed. 

The  problem  of  approximating  q  is  of  the  same  nature  here 

as  in  the  original  formulations  of  the  integral  equation  methods. 

But  there  one  has  sor  a  simplifying  features.  The  boundary  conditions 

for  the  profile  are  prescribed  along  the  x  axis.  In  the  case  of 

a  symmetric  flow  field  the  solution  for  the  linearized  problem  is 

directly  expressed  by  quadratures.  Unsymmetric  flows  are  more 

complicated  but  at  least  the  linearized  flow  field  can  be  computed 

once  and  for  all.  In  the  present  case  the  inner  contour  of  the 

distant  field  is  likely  to  be  a  rectangle.  Along  this  boundary 

neither  the  normal  derivat.ive  nor  the  potential  is  known.  By 

the  computation  carried  out  in  the  inner  field  one  obtains  at  the 

inner  boundary  of  the  distant  field  only  a  relation  between  the 

two  quantities  (or  other  related  data,  for  instance,  the  velocity 

components  u  and  v) ,  and  this  relation  changes  from  iteration  to 

iteration.  The  relation  so  obtained  must  be  combined  with  the 

far  field  conditions.  Then  one  obtains  o  and  4  or  u  and  v  at 

Tn 

the  boundary.  These  are  the  quantities  on  which  the  estimation 
of  q  in  the  distant  far  field  must  be  based. 

At  this  stage  the  experience  gained  with  the  classical 
integral  equation  method  may  be  helpful.  To  give  an  idea  of 
what  one  might  do,  1  describe  here  a  procedure  which  has  much  in 
common  with  Nixon's  extended  integral  equation  method.  For 
simplicity  we  assume  that  the  inner  contour  of  the  distant  field 
is  given  by  two  lines  y  =  constant  at  some  distance  above  and 
below  the  profile  (Fig.  12).  The  computed  part  of  the  flow  field 
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would  then  be  a  strip  in  the  x,y  plane.  In  practice  one  will, 
of  course,  cut  off  the  strip  by  lines  x  =  constant;  then  one 
needs  additional  conditions  for  these  lines.  We  restrict  ourselves 
to  the  simpler  problem  and  derive  expression  for  the  upper  part  of 
the  distant  field. 

Let  y^  be  the  value  of  y  for  the  inner  boundary  of  the 
distant  field.  Assume  that  u  is  known  at  y  =  y^  and  at  another 
line  y  =  y2;  Y 2  ;>  Y]_*  Along  lines  x  =  constant,  u  is  then 
approximated  by 


u  (x,y)  = 


2 T! 


2  2 
(y  +x  ) 


a  (x) 


b  (x) 
3 
Y 


The  first  term  is  due  to  the  circulation,  the  values  of  a(x)  and 
b(x)  must  be  determined  in  such  a  manner  that  u  assumes  the 
desired  values  at  the  lines  y  =  y^  and  y  =  y2 •  It  is  consistent 
with  the  character  of  the  approximation  if  we  simplify  the  last 
equation  to 


u(x,y)  =  ^7  y  1  +  a(x)y  2  +  b(x)y  3 

of  course,  with  modified  functions  a(x)  and  b(x).  One  obtains 
u(x,y)  =  F/  ( 2 7T y )  (41) 

+  (yl"y2)'1  f[k(x'Yi>  "  C  P/ (2Tty1)  )  )yi3-(u(x,y2)-  (  V/ (2^)  )  )  y.,3]  y 
+  [- (u  (x,yx)  -  (  V  ( 2 Tty  1 )  )  )  y3y2  +  ( u  (x  ,y2 )  -  ( 1/ ( 2ny2 )  )  )  y^y 3  ]  y-3 } 

as  can  easily  be  verified.  For  other  contours  one  might  think  of 
approximations  of  a  similar  character. 

Assume  now  that  we  carry  out  a  procedure  based  on  Eq.  (34). 
There  one  carries  out  an  iteration  for  the  computed  part  of  the 
flow  field  in  which  the  function  f(s)  changes  in  each  iteration 
step.  We  want  to  determine  an  approximation  to  q(x,y)  based  on 
Eq.  (41)  which  pertains  to  a  fixed  function  f(s).  This  is  done 


iteratively.  Assume  that  one  has  some  approximation  to  q.  One 
then  uses  Eq .  (34)  to  evaluate  u  for  y  =  y  ^  and  y ^ •  One  has 
specifically: 

u  ( x , y )  =  (2  0  1 (  -  f  f(u)yr(f(o)  -  x,  y  -y)do 

(a: 

-fS  q  ( ■■  ,  n)  (t.-x,  n-y)  df;dn  +  ir(x,y)  '  sinu(x,y)] 


These  values  of  u  are  then  substituted  into  Eq .  (41),  and  with 

this  equation  values  of  q  throughout  the  distant  field  are 
defined.  With  these  values  the  computation  is  repeated.  Of 
course,  with  a  more  complicated  inner  boundary  of  the  distant 
field,  the  formulae  for  computing  u  at  the  pivotal  lines  of  the 
distant  field  will  become  somewhat  more  complicated  and  one  will 
in  some  regions  use,  instead  of  Eq.  (41) ,  somewhat  different 
interpolation  formulae. 

It  is  likely  that  only  very  few  iteration  steps  will  be 
needed.  The  values  of  q  so  obtained  are  then  used  in  the 
iteration  process  for  the  computed  part  of  the  flow  field  des¬ 
cribed  in  the  preceding  section.  The  key  formula  is  Eq.  (38). 

The  form  of  the  far  field  conditions  just  discussed  is 
probably  best  for  steady  flows  because  they  yield  local  boundary 
conditions  for  (or  4>t)  at  the  outer  boundary  of  the  computed 
flow  field. 

Because  of  their  global  character  the  far  field  conditions 
Eq.  (31)  can  be  used  only  for  linear  problems;  for  the  number 
of  simultaneous  equations  which  one  must  solve  is  too  large. 

Such  linear  problems  arise  if  one  tries  to  improve  a  sufficiently 
close  approximation  to  the  solution  by  a  Newton-Raphson  procedure, 
or  if  one  computes  small  periodic  disturbances  to  a  steady  flow 
field.  In  the  latter  case  one  may  be  forced  to  use  these  global 
conditions  if  the  frequency  is  too  high  or  the  Mach  number  is 
too  close  to  1.  (The  question  of  how  important  such  cases  are 
for  technical  applications  is  left  open.)  We  sketch  the  procedure 
applicable  for  these  global  boundary  conditions. 
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One  chooses,  as  solution  strategy  for  the  inner  region,  a 

procedure  which  yields,  actually  as  an  intermediate  step,  relations 

between  the  potential  and  its  normal  derivative  at  its  outer 

boundary.  These  relations  are  combined  with  the  far  field 

conditions  in  their  global  form,  Eq.  (31)  .  At  this  stage  one 

must  use  some  approximation  for  q,  for  instance  q  =  0.  One 

then  obtains  the  values  of  the  potential  and  its  normal  derivative 

or  rather  the  velocity  components  u  and  v  at  the  outer  boundary 

of  the  computed  field  (which  is,  of  course,  the  inner  boundary 

of  the  distant  field).  One  then  uses  Eq .  (A17)  (with  approximate 
2 

values  of  u  in  the  field  to  compute  u  for  y  =  y  ,  Eq.  (A17)). 

2nu(x,y  )  +  /  u(f,,y,)ip  U>x,y  -y  )df 

Z  'Y'l  -L  r)  l  z 

Ci4,2 

2 

+  /  v  ( f, ,  y , )  [  v  .  ( f,-x ,  y ,  -y  9 )  dC  -(n/2)u  (x,y_) 

-  (1/2)  //  u2(f',n)  /rr  (  r-x  ,  r]-y)  dfdg  . 

ft 

The  values  so  obtained  u(x,y2>  together  with  the  values  of 
u(x,y^)  with  which  this  step  starts  are  then  substituted  into 
Eq.  (41)  and  the  expression  for  u  so  obtained  is  used  in  the 
global  form  of  the  far  field  conditions,  Eq.  (31).  Because  of 
the  changes  in  q  which  occur  in  this  procedure,  one  will  now 
obtain  new  values  of  u  and  v  at  the  outer  boundary  of  the 
computed  field  and  with  these  values  the  approximation  to  the 
values  of  u  in  the  distant  field  is  improved. 

We  mentioned  above  that  the  procedure  will  be  applied  either 
to  linear  or  linearized  problems.  Accordingly,  there  will  be 
a  linearization  also  in  the  far  field  conditions.  Ic  will 
depend  upon  the  special  circumstances  whether  a  direct  or  an 
iterative  approach  is  then  preferable. 
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SECTION  IX 


CONCLUDING  REMARKS 

The  report  explores  the  integral  equation  method  from  a 
theoretical  point  Oi  /iew  in  some  depth,  but  without  losing  sight 
of  its  practical  purpose.  It  is  shown  that  the  integral  equation 
method  can  be  regarded  as  a  method  of  weighted  residuals  with 
rather  unconventional  weight  functions. 

The  resulting  equations  contain  terms  obtained  from  linear 
theory  which  are  solely  determined  by  data  at  the  boundary  of 
the  region  under  consideration  and  further  nonlinear  terms  which 
depend  upon  flow  field  data.  These  latter  terms  are  of  smaller 
significance  and  can  therefore  be  evaluated  with  reduced  accuracy. 

In  the  original  form  of  the  integral  equation  method  these  field 
terms  are  approximated  by  means  of  boundary  data.  This  method  is 
particularly  valuable  for  purely  subsonic  flows. 

In  the  transonic  region,  nonlinear  flow  field  data  from  the 
interior  gain  importance.  To  evaluate  these  data  Nixon  uses 
again  the  formulation  of  the  integral  equation  method.  Since  the 
primary  goal  is  to  make  the  residuals  in  the  partial  differential 
equation  small  in  the  region  close  to  the  profile,  it  appears 
preferable  to  use  these  formulations  which  use  local  weights  rather 
than  the  global  weights  of  the  integral  equation  method.  Accordingly, 
the  present  paper  advocates  hybrid  methods  in  which  the  integral 
equation  formulation  is  used  at  points  where  the  nonlinear  terms 
are  of  small  importance  while  finite  difference  or  finite  element 
methods  are  used  in  the  region  close  to  the  profile.  Such  a 
procedure  makes  it  possible  to  formulate  the  shock  conditions  and 
the  conditions  at  the  sonic  line  in  a  local  form,  which  is  closer 
to  the  nature  of  the  problem  at  hand.  The  form  in  which  these 
conditions  are  expressed  by  the  integral  equation  method  is  very 
ingenious,  it  is  true,  but  still  somewhat  awkward. 

In  an  extreme  development  one  uses  the  idea  of  the 
integral  equation  method  only  in  the  distant  field  in  order  to 
derive  boundary  conditions  at  the  boundary  between  the  distant 
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field  and  the  computed  part  of  the  flow  field.  The  far  field 
conditions  found  in  the  literature  assume  that  in  the  distant 
field  the  linearized  potential  equation  holds.  The  method 
proposed  here  borrows  from  the  integral  equation  method  the 
idea,  that  the  information  about  the  distant  field  that  is  needed 
to  take  nonlinear  terms  into  account,  can  be  inferred  from  the 
data  for  potential  and  its  normal  derivative  at  the  boundary 
between  computed  flow  field  and  distant  field.  One  can  hope  that 
such  improved  boundary  conditions  allow  one  to  reduce  the  size 
of  the  computed  part  of  the  flow  field  without  undue  loss  of 
accuracy . 

The  development  has  been  carried  out  for  two-dimensional 
problems.  But  extensions  to  the  three-dimensional  problem  are 
not  difficult.  Such  extensions  have  already  been  made  by  Nixon, 
but  of  course  not  for  the  hybrid  form  of  the  method  and  the  far 
field  conditions  advocated  in  the  present  report. 


f 
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APPENDIX  I 


DERIVATION  OF  THE  SHOCK  CONDITIONS  FROM  THE 
POTENTIAL  EQUATION  BY  MEANS  OF  THE 
CONCEPT  OF  WEAK  EQUALITY 


The  starting  point  is  the  potential  equation  (5) 

fx  <*x  -  (*x2''2>>  +  3^  <y  -  <A1> 

Let  w(x,y)  be  a  weight  function  of  finite  support  (this  means  that 
w  is  different  from  zero  in  a  finite  region  only) .  It  is  assumed 
that  all  derivatives  which  occur  in  the  analysis  exist  and  that  w 
and  its  first  derivative  vanish  at  the  boundary  of  the  support  region. 
The  postulate  that  Eq.  (Al)  is  satisfied  in  the  weak  sense  is  ex¬ 
pressed  by 

//  w(x,y)  [  -  Ux2/2))  +  ^  (4>y)]  dxdy  =  0 

for  every  function  w(x,y).  Carrying  out  the  usual  treatment  of 
derivatives  for  generalized  functions,  one  obtains  by  means  of  in¬ 
tegrations  by  part 

‘  //I*x  '  <*x2/2»  U  +  *y  fjbdxdy  =  0. 

The  integral  over  the  boundary  of  U  vanishes  because  of  the  defini¬ 
tion  of  w.  The  integrand  is  bounded  even  at  the  shock.  The  expres¬ 
sion  therefore  remains  unchanged  if  one  excludes  the  shock  by  a  strip 
of  width  2r  and  then  lets  e  tend  to  zero.  Let  the  region  so  modified 
be  denoted  by  £2’  (Fig.  1). 

In  the  region  q'  the  partial  differential  equation  is  satisfied 
pointwise.  Now  we  carry  out  an  integration  by  parts  in  the  last 
equation.  The  boundary  now  includes  portions  and  which  ex¬ 

tend  upstream  and  downstream  of  the  shock  in  the  distance  e.  Let 


the  shock  be  given  by  x  =  xs(y).  Then  one  has  for  9ftj_ 
x  =  xs(y)  -  £ 
and  for  9ft  .j 

x  =  xg (y)  +  e  . 

The  double  integral  which  arises  by  the  integration  by  parts  vanishes 
because  <J>  satisfies  Eq.  (Al)  point  wise  in  ft '  .  The  contribution  of 
the  outer  contour  of  ft'  vanishes  because  of  the  definition  of  w. 

One  therefore  obtains 

/  w{  [<J-X  -  (4>x2/2)3  dy  -  <(>y  dx]  =  0 
9ft^+9ft2 

The  direction  of  integration  is  to  be  chosen  so  that  the  region  of 
integration  lies  to  the  left.  This  means  that  one  proceeds  along 
9ft^  and  3ft2  in  opposite  directions.  One  has,  because  of  the  defini¬ 
tion  of  w 

lim  (w ( x  ( y )  -  e,y)  =  lim  w(x  (y)  +  e,y)  =  w(x  ,y) 

£-*■0  S  e->0  S 

Using  the  notation  of  Section  hi  one  therefore  obtains 

/  w(x  (y),y)  {[4.  -  (<J>  2/2 )  ] +  dy  -  [<f>  }+  dx}  =  0 

9ftx  S  X  x  y 

In  order  for  this  equation  to  be  correct  for  any  choice  of  w,  one 
must  have 

?  +  dx 

[»x  -  [V-  -w  =  0  • 
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APPENDIX  II 
THE  EXPRESSION  Q(x,y) 


We  have  defined  in  Eq.  (21) 

Q(x, y)  =  If  [yl  (<Pr-(4>£2/ 2))  +  ($n)  ]i|»(5-x,n-y)dCdn 


At  the  moment  ft  is  some  region  of  the  f>fn  plane  in  which  ip  satisfies 
the  Laplace  equation.  Its  boundary  is  denoted  by  3ft.  The  region  may 
include  shocks.  One  obtains  in  a  familiar  manner  by  carrying  out 
integrations  by  part 

Q(x,y)  =  /  iM£-x,o-y)  [<J>r(£,n)dn  -  <J>  U,n)d£;] 

3ft  t,  n 


-  j  41  (C,n)  [>Jv  (£-x,n-y)dn  -  i|>_U-x,n-y)d£;] 
3  ft  n 


-  (1/2)  fj  <m-x,n-y)  (^2(t,n)dWn  (A2) 


A  further  double  integral  with  integrand  <J>  (iJj^  +  ij>  )  vanishes  because 
ip  satisfies  the  Laplace  equation.  The  integration  around  the  contour 
3ft  is  performed  in  the  direction  for  which  the  interior  of  ft  lies  to 
the  left.  One  has 

-  (1/2  ff  iM£-x,n-y)  -Jy  (<Pr2)  d^dn  =  -  \  j  i|»(£-x,n-y)  <J>r2dri 
ft  ^  3D  ^ 

2  (A3) 

+  (1/2  //  <p  (i-x,n-y)<Pr  (F.,n)dCdn 
ft  ^  ^ 


The  form  of  the  right  hand  side  in  Eq .  (A3)  may  sometimes  be  pre¬ 
ferable  because  it  eliminates  one  differentiation  in  the  unknown  func¬ 


tion  <P“ 

2, 


It  thus  becomes  applicable,  even  if, in  some  approximation, 

The  integrands  of  the  double 


3  ((j)._  *■)  /  3^  contains  delta  functions, 
integrals  are  singular  at  the  point  (£,n)  =  (x,y) 


43 


We  choose  for  q  a  region  of  the  £,n  plane  punctured  at  the 
point  ( £ ,  n )  =  (x,y)  (Fig.  2).  Its  boundary  is  3£2.  Part  of  3Q 
at  a  great  distance,  to  be  denoted  by  3 S2 ^  is  given  by  a  circle 
with  a  radious  R.  Subsequently,  the  limiting  process  R-*-^  will 
be  made.  In  the  familiar  forms  of  the  integral  equation  method 
the  profile,  extending  from  £  =  0  to  £  =  L,  is  replaced  by  a  cut 
along  the  r-axis.  To  prepare  for  a  modification,  we  allow  for 
the  possibility  that  the  part  of  3P  which  is  usually  given  by  the 
two  borders  of  this  cut  lies  at  some  distance  from  the  profile. 

This  part  of  3P  is  denoted  by  3P„.  The  wake  is  excluded  by 
a  cut  along  the  £-axis.  The  two  borders  of  this  cut  form  the 
boundary  They  are  joined  to  either  at  the  trailing  edge 

£  =  L  or  at  some  point  further  downstream  £  =  .  The  junction 

with  3£>  lies  at  £  =  R.  The  small  circle,  by  which  the  point 

( £  ,  n )  =  (x,y)  is  excluded  is  denoted  by  3P  . 

In  the  presence  of  a  circulation  the  term  of  the  potential 
which  are  dominant  at  a  large  distance  have  the  form 

<(>  =  const  -r0/2ir 

where  F  is  the  circulation  and 

Q(£,n)  =  arc  tg  (ri/£)  ;  0  £  2it. 

The  choice  of  the  constant  in  the  above  equation  has  no  effect  on 
the  flow  field,  because  the  flow  differential  equation  contains  only 
derivatives  of  <j> .  Accordingly,  we  can  choose  (for  R  large) 

<j>  =F  (  (1/2)  -  ( 0/2tt )  )  .  (A4) 

As  R  tends  to  infinity,  <{>  and  tend  to  zero  strongly  enough 

that  for  the  boundary  3p^  the  first  integral  on  the  right  of 

Eq.  (A2)  vanishes.  Regarding  the  second  integral  one  observes  that 

(J< £- dig  -  <Jjfid£  =  (dif;/dn)do 

where  d/dn  denotes  the  derivative  in  the  direction  of  the  outer 
normal  and  do  is  the  arc  length  of  3fi.  For  R  sufficiently  large 


and  (::,y>  bounded,  one  has 


ip  -  log  R 


do  =  RdO. 


One  then  obtains 

-  2  TT 

-  /  0  (dtp/dn)  ds  =  T  (0/2)  -  (0  /4tt)  J  =  0. 

3ft.  0 


The  contribution  of  the  contour  3ft^  vanishes,  if  one  chooses  the 
constant  in  <p ,  which  is  arbitrary,  according  to  Ec.  (A4). 

No  simplifications  arise  in  the  integral  for  a  general  choice 
of  the  portion  of  the  contour  3Q2-  In  cases  where  3ft»  is  represented 
by  the  two  borders  of  a  cut  along  the  x  axis  (from  f=0  to  £=L)  one 
can  treat  3ft2  anc^  ^3  together.  On  the  lower  side  of  the  cut  one 
moves  in  the  negative,  on  the  upper  side  in  the  positive  £,  direction. 
In  this  case  dg=0  and  one  obtains 


/  ip(£-x,n-yH  (£,n)d£  =» 

9ft2+3ft3  1 

L 

-  (  ip  ( £-x , -y )  { <p  (  5  > +0)  _  ^  (  £ ,  _0)  }d£ 

0  n 

Here  we  have  written 

lim  <J>n(C,n)  =  <Pn»£,+0),  lim  <J>nU,n)  =  <J>n  ( 5 ,  —0 ) 

n-’-o  n-*-o 

ri>0  n<0 

The  wake  does  not  appear  in  the  integral  because  there  <j>  is  con¬ 
tinuous.  At  the  upper  and  lower  side  of  the  profile,  <J>  is  given  by 
the  boundary  conditions.  For  a  symmetric  profile  one  has 

4>n  (5,-0)  =  -<t>n  ( 5,  +0)  . 
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If  3^2  +  3^3  given  by  the  two  borders  of  a  cut  along  the  axis, 
then  one  obtains  for  the  second  integral  in  Eq .  (A2) 


/  'P(C,0)<i'  (Z-x,n-y)dFJ  =  I  [<{  ( £ , +0)  -  <p  (  £ ,  -0)  ]  >}>  (  £-x,  -y)  d£ 

3fi2+3H3  1  0  n 


At  the  profile  the  difference  tp(£,+0)  -  <M£,-0)  is  not  known.  For  a 
symmetric  flow  field  it  vanishes  along  3SI_  +  In  t^ie  wa^e 

this  difference  is  given  by  the  circulation.  Substituting,  also 
the  expression  for  ip ,  Eq .  (  13),  one  obtains 


+  r  /  ip  (£-x,-y)d£  =  -r  n  (x-r  ,y)  J 
L1  L1 


]'  [0  (x-L^  ,y)  -it] 


where 


0(x-£,y)  =  arc  tg(y/(x-f,), 

0  <  0 (x-£ ,y)  <  2  tt  . 

The  contour  integral  around  the  singular  point  is  treated  in 
a  familiar  way.  Let 

£-x  =  r  cos  0 
n-y  =  r  sin  0 

The  first  integral  over  in  Eq .  (A2)  vanishes,  because  the  singu¬ 

larity  ip  =  log  r  is  too  weak.  For  the  evaluation  of  the  second 
integral  we  note  that 

ib  ^  =  r  ^  cos  ,  ib  =  r  ^  sin  0 

C  n 

df  =  -r  sinO  ip  =  r  cos  0 

T  n 


fore 


The  integration  proceeds  in  the  negative  0  direction.  There- 


lim  (-  /  <)>(£, n)  [ipr  (f-x,n-y)dn  -  ipr  (C-x, n-y)  3  )  =  2n<p(x,y) 
r  ►o  an ,  s  n 
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One  thus  obtains 


Q(x,y)  =  2M'  (x,y)  + 


-x, n-y)  [  i  .  d n  -  d' 

f  r 


2 


-  /  v  (/ ,  n)  [or  (t-x,  n-y)  dn  -  t  ( '-x,  n-y) d;. ; 

n  '3  '! 


(A5) 


- !’  ( tt-0  (x-l,  , y )  -(1/2)//  y  (/-x,  n-y)  rryr( 'P  ~  ( / ,  n )  d/,dn 

si  ‘  f  - 


Eq.  (22)  is  obtained  by  specializing  rh';-  to  a  cut  along  the  x 
axis . 

A  special  limiting  process  familiar  from  potential  theory 

is  needed  if  the  point  x,  y  approaches  a  point  x  ,  yc  of  the 

contour  along  the  inner  normal.  One  then  treats  a  vicinity  of 

the  point  xc»yc  separately.  Let  be  a  small  portion  of  the 

contour  SO.  which  contains  the  point  x  ,y  ,  and  let  = 

c  c 

if  For  (x,y)  =  (xc,yc)  the  integrand  in 


f  4>(trn)  [/c(i-x  , n-y_)dn  -/,(:-  x  ,  n-yjds 

3H"  c  c  c  c 

remains  bounded  provided  that  3/  has  finite  curvature  at  (x  y  ) . 
therefore  can  immediately  make  the  limiting  process  in  which  the 
excluded  portion  tends  to  zero.  Discussions  familiar  from 
potential  theory  give 


lim  /  4>(,r.,n)  [if- (t-x, n-y) d/  -p  (f,-x,n-y)dn  =  "<f>(x  ,y  ) 

(x,y)-(xc,yc)  DR*  ‘  n  c  c 

Including  also  cases  where  the  point  (x,y)  lies  inside  of  ^ 
we  rewrite  Eq.  (A5)  in  the  form 


Q(x,y)  =ot  tt  <fi  (x,y)  +  /  ;j<  ( f,-x ,  n-y  )[<f>r(  t ,  n)  dp  -<t>n  ( /  r  n)  d;j 

a‘;2  •  ” 

-  / ( Z , n )  [iir  (f.-x, n-y) drt  -ip  (£-x,n-y)d/]  (A6) 

mi  -  n 

2 

-  r  (  n  — 0  (x-L  ,  y ))-  ( 1/2 )  /.A,;  (f-x,n-y)~ (4>J  (/,n)  )d\dn 
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We 


<  =  1  for  (x,y)  on  0^ 
i  =  0  for  (x,y)  outside  n. 

One  has  (:n  the  second  integral  over  ^  ) 

lim  {>.  ( •'  (a)  -  x  ( s )  ,n  (o)-y(s)  )^(o) 

X  (s)  a° 

y(s)  -r.  (o) 

-tj>,  ( '  .  (o  )  -x  (s)  ,  r,  (o)  -V(s))  ■  =  (1/2)  R(o)_1 

where  R(c)  is  the  radius  of  curvature  at  (o),n(o),  positive  if 
is  convex. 

The  expression  Q(x,y)  is  continuous  along  a  line  y  =  const  as 
the  point  (x,y)  passes  through  a  shock.  This  is  seen  in  the  follow¬ 
ing  manner.  We  have  postulated  continuity  of  h  In  the  integrals 
over  v  and  its  derivatives  are  continuous,  because  the  shock  is 

not  part  of  and  therefore  (x,y)  /  (f,n).  A  discontinuity  can 

therefore  be  introduced  only  by  the  double  integral.  It  would  be 

2 

caused  by  the  jump  of  oil  at  the  point  (C,n)  =  (x,y)  if  the  point 
(x,y)  passes  through  a  shock.  For  a  discussion  is  suffices  if  one 
considers  the  critical  integral  for  the  vicinity  of  some  point  of 
the  shock  x  =  xs  (y)  with  <J>f  =  1  upstream  and  =  0  downstream  of 
the  shock.  We  enclose  this  point  by  a  rectangular  contour  (Fig.  3) . 
Distance  points  of  fi  will  not  give  a  discontinuity  in  Q,  therefore 
it  sufficies  if  one  evaluates  the  integral 

( 1/2 ) f  ( c~x  > n-y ) d ' dq 

once  for  (x,y)  =  (x  -  e,y)  and  a  second  time  for  (x,y)  =  (x  h  e,y) 

s  s 

and  examine  the  difference  in  the  limit  e  -  0.  The  integration 
is  to  be  carried  out  for  the  shaded  area  of  Fig.  3.  One  obtains, 
by  carrying  out  the  integration  with  respect  to  £ 


n2  n2 

(1/2)  [-  f  ip  ( f.1~xs+e ,  n-v)  dn  +  /  i Hf,  (n)-x  +c,n-y)dn 

nl  ni 

in  the  second  case  with  e  replaced  by  -e.  The  first  integral 

which  comes  from  the  left  boundary  of  the  rectangular  box,  will 

not  cause  a  singularity  as  the  sign  of  c  changes.  In  the  second 

integral  we  place  the  origin  of  F, ' ,  n*  system  into  the  point 

( x  , y )  and  approximate  F  '  (n)  by  e'  =  -  tgBn’,  where  3  is  the 
b  s  s 

angle  of  the  normal  to  the  shock  at  the  poJ.nt  in  question  with 
the  x  axis.  Finally  we  substitute  p ;  one  then  obtains 

1 

n2 

-  (1/4)  /  log  (( —  n  '  tgB+e )  2  +  n' 2)  dn! 
nl1 

This  integral  can  be  evaluated  in  closed  form.  The  resulting 
expression  is  continuous  as  c  passes  through  zero. 
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APPENDIX  III 
DERIVATIVES  OF  Q(x,y) 


After  replacing  the  last  term  by  means  of  Eq.  (A3)  one  ob¬ 
tains  from  Eq.  (A5)  by  differentiation  with  respect  to  x 

Qv(x,y)  =  2tt6  ( x , y)  +  /  -\p f(£-xf  n-y)  [<j>,dn-<)>  d£  ] 


-  /  <P(f,,n)  l-4>,  r(C-x,n-y)dn  +  iK  (£-x,n-y)d£) 

9fi2  *>■•> 

oo 

-  r  /  ^^n(C-x,-y)dC 

L1 

-  (1/2)  /  -  ip  (C-x,n-y)  ip,2  (t,n)dn 

dn2 

+  (1/2)  ^  //  <^(£-x,n~y)c^2(C,n)d£dn  (A7) 


In  the  expression  with  the  factor  Y  we  have  substituted  the 
form  in  which  it  originally  appeared  in  Appendix  II.  One  can  carry 
out  the  integration  with  respect  to  £  and  then  obtains 

oo 

-r  /  *  (£-x,-y)d*;  =  r  ^(Lj-Xf-y)  (A8) 

L1 

In  the  second  integral  over  3ft2  we  use  the  relation 
For  the  differentiation  along  a  contour  given  by 

5  -  «c(n) 

an  operator 
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D  ^  d£  » 

c  _  9  c  9 

dn  9n  +  dn  9^ 


is  introduced.  One  then  obtains  for  the  second  integral  over  9fi„ 


/  <P(£,n)  [-^j.^(C-x,n-y)dn  +  (£-x,n-y)d£  ] 


-  /  (it»n)dn  =  -  4>U#n)<l>n 


|  +  /  %  Hn(<M5'n)dn 

limits  dn2 


=  -<t>U»n)<f'  I  +  /  i|»  (<t>  (C,n)  dn  +  <j>,  (S,n)dC) 

'  .  1  . .  9ft0  1  1  ^ 

limits  2 


The  contour  9^  starts  and  ends  at  the  point  £  =  ,  n  =  0  where 
the  two  borders  of  the  wake  contour  90^  start  or  end.  There  the 
terms  for  the  limits  are  canceled  by  the  expression  Eq .  (A8).* 

We  write 


4>X  =  U,  <Py  =  V 


The  integrals  over  '<> P ^  i-n  Fq*  (A7)  then  contract  to 


(A10) 


{  -(u(C»n)  -  (1/2) u  (f,, n))  <j>r  (C-x,n-y)dn  +  u(C»n)^  U-x,n-y)dE; 


+  /  vU,n)  [4v  ( C-x,  n-y) dC  +  ^  (f>x,n-y)dn] 

9^2  ’  1 


*This  simplification  occurs  only  for  steady  plane  flows. 


(All) 


i 


In  evaluating  the  derivative  of  the  double  integral  in 
Eq.  ( A7 ) ,  some  precautions  are  necessary  because  of  the  singularity 
of  the  integrand.  In  the  original  definition  of  the  region  2 
the  singular  point  (£,n)  =  (x,y)  is  excluded  by  a  small  circle 
of  radius  e  which  moves  as  the  point  (x,y)  moves.  This  must  be 
taken  into  account  in  the  differentiation  with  respect  to  x  or  y 
(up  to  now  the  exclusion  of  the  point  has  not  mattered  in  the 
evaluation  of  the  double  integral) .  We  divide  the  region  Q  into 
two  regions  by  surrounding  the  point  (x,y)  by  a  box  of  finite 
dimensions  (Fig.  4) .  The  boundaries  of  the  box  remain  fixed 
while  the  differentiation  with  respect  to  x  or  y  is  carried  out. 

In  the  outside  region  the  differentiation  can  immediately  be 

carried  out  under  the  integral  sign: 


(A12) 

//i|JC-x,r)-y)  u2  ( c / n)  d£dn=-(l/2)  //  ip  ,  (C-x,  n-y  (u2  (C  ,n)  dCdn 

n,  n,  ^ 

X  x 


For  the  evaluations  within  the  box  (region  ft2)  we  set 


?1  = 

nl  = 


<v~X 


n-y 

-2 

•1 


(A13) 


r  =  (C7  +  n 


2.1/2 

r 


The  boundaries  of  the  region  are  given  by 


x^ , and 

x2. 

xl 

<  X  < 

X2 

ylfand 

n=  y2, 

yl 

A 

A 

y2 

and  the  small  circle  with  radius  c  around  the  point  (x,y) .Let 
temporarily 

I  =  (1/2)^|  /  /  (C-x,n-y)u2(C,n)dCdn 

n=y1 

r  >e 


then,  after  substitution  of  Eqs.  (A13) 


52 


x2-x 


y,-y 


=  (1/2)  ^  Him  / 

e->0  £,=x.-x 

±  X 


/  (^n-^u  t^+x, nj+yjd^dri  ] 

ryry  ' 


In  the  £.  n,  system  the  small  circle  around  the  singular  pcint 
11  2 
is  fixed,  even  if  x  or  y  changes.  The  function  u  depends  upon 

two  arguments  (originally  £  end  n) .  For  clarity,  we  denote 

differentiation  with  respect  to  its  first  or  second  argument, 

respectively,  by  subscripts  1  and  2.  Now  we  carry  out  the 

differentiation  with  respect  to  x  in  the  last  expression 


y-,-y 


y  ~y 

I_  (i/2)T  J;  'J'f ;(x1-x,  H-^u2  (x1,n1+y)dn1-  J  lpjx2~x/  )  u2  (x2  ,  Hj+y )  dr^ 

L  n.  =v _  —  v  r. .  =v.-v 


n1=y1~y 


X2"X 


ni=yry 


y0-y 


+  lirn  /  /  i|/  (£  ,  H-.)  u2  (C-.+X,  ru+y)  dC.dn, 

C1=x1-x  ri1=y1-y 


Returning  to  the  original  variables  £  and  n  one  obtains 

y-,  y2 

“22 
I  =(1/2)  [  /  ip  (x.-x,n-y)u  (Xj^/TiJdn  -  f  p  (x2~x,r)~y)  u  (x2,n)dn 

n=yi  J  n=y-,  ^ 

1  (A14) 

x2  /  2 

+  lim  /  /  (£-x,n-y)  O  (u2  (£,n)/9£)d£dn  ] 

e+0  £=x1  n=y-|  ** 

In  the  double  integral  an  integration  by  parts  with  respect  to  £ 

1  2  2  2 
is  carried  out.  Let  9fi2  be  the  circle  (£-x)  +  (n~y)  =  £  . 

One  obtains 


c  2  c  2 
ff=  -  f  ipr  (Xj^-x,  n-y) u  (Xj^nJdn  +  /  <p  (x2-x,n-y)u  (x2,n)dn 

n=yx  ^  n=y1 

(A15) 

lim  /  (|;  (£-x,n-y)u2(£,n)dn  -  lim  //  f(£-x,n-y)u2(£,n)drjdn 

£->•0  an,  •  r-o  u  '*'• 

rN0  2 


.vJro 


kW-  v.; 


The  first  two  integrals  in  Eq.  (A15)  cancel  the  first  two  integrals 
in  Eq.  (A14) .  We  introduce  temporarily 


r  - x  =  r  cos  0 
n-y  =  r  sin  0 

I 

Then  along  dQ^ 

do  =  e  cos  0  do 


Moreover 

ip^  -  c  cos  e 

The  integration  around  proceeds  in  the  negative  B  direction. 

One  obtains 

2  2 

lim  /  iji  (f,-x,n-y)u  (fJ/n)dn  =  -  nu  (x,y) 
e-*0  3Q' 

The  double  integrals  in  Eqs.  (A15)  and  (A12)  can  be  combined.  One 
finally  finds 

(1/2)  ~  [  ff  <//,  (c-x,  n-y)  u2  (i' ,  n)  d^dr)  = 

n  •  (A16) 

-u/2u2  (x,y)  -  (1/2)  If  i>  (C-x,n-y)u2  (C,n)dCdn 

p  ^  * 

Substituting  these  results  into  Eq.  (A7)  one  obtains 

Qx  ( x , y )  =  2  ttu  (x,y)  +  (A17) 

/  -  (u(£ ,n) - (1/2) u2  (£, n)  ] <Pr  (S-x, n-y) dn  +  u(£#n)<l>  (£;-x,n-y)d£ 

+  /  v(^,n)  [>pr  (Z-x,  n-y)  dr,+p  (f,-x,n-y)dn  ) 

oa2  5 

-  (rr/2)  u2  (x,y)  -  (1/2)  // u2  (C.n)  U-x,n-y)  d£dn 


Specializing  this  equation  to  the  usual  case,  where  is  given 
by  the  two  borders  of  the  cut  made  along  the  profile  (from  x  =  0 
to  x  =  L) ,  one  obtains: 


L 

Q  (x,y)  =  2ttu  (x,y)  +  f  [  (u  (£ , +0 ) -u (£  ,-0)  ]  (£-x, -y)  d£ 

x  0 

L 

+  /  [v (£  ,+0)  -v( £ , -0 ]  <Jj  r  ( £-x,  -y)  d£  _ 

0 

-  (tt/2)  u2  (x,y)  -  (1/2)  // u2  (£,  n)  (£-x,n-y)  d£dn 


A  corresponding  procedure  can  be  carried  out  for  Q  . 
For  the  specialized  region  U  one  obtains 

L 

Qv(x,y)  =  2ttv (x ,y)  +  /  [v (£,  +0)  -v ( £ , -0)  ]  ( £-y , -y ) d£ 

Y  0  1 


-  /  (u ( £  ,+0)  -u  (£»  -0)  (£-x,  -y)d£ 

0  ^ 

~ (1/2)  //  u2 (£, n)i|vT1  (£-x,n-y)d£dn 

Q  K n 


(A18) 


( A19 ) 


APPENDIX  IV 


DERIVATION  OF  THE  EXPRESSION  Q  (x,y)  BY  APPLYING 
THE  WEIGHT  FUNCTION  -pr  TO  THE* POTENTIAL  EQUATION 

The  weight  function  -ip^  has  a  singularity  at  the  point 
(t',n)  =  (x,y).  The  question  arises  to  what  extent  this  fact 
imposes  a  restriction  on  the  form  of  the  residuals  to  which  this 
weight  function  can  be  applied.  Introducing,  as  before 

£-x  =  rcosO,  n-y  =  rsinB 


one  has 


=  -r1  cos  6,  =  -r1  sin6 

We  postulate  that  the  residuals  to  which  this  weight  function  is 
applied  be  continuous  and  that  its  first  derivative  exist  at  the 
point  (x,y) .  Such  a  residual  can  then  be  written  as 


R(£,n)  =  Rq  +  R^.  rcosQ  +  R^  r  sin0  +  terms  of  higher  order 
Here  R^,  R„  and  R^  are  constants.  One  then  considers  expressions 
-  //  ( Rq  +  R^  r  cos0  +  R^rsinO)  ip ^  d£  dn 


r>e 


-  "  //( Rn  +  R,  r  cos0  +  R  rsin0)r  ^ cos0  r  dr  d0 
u  i  n 

r>t 

0<_0  <2v 

2  tt  2  2 

_  _  j  [r  Rq  cosO  +  (r  /2)  ( R^  cos  0  +  sinO  cos0)] 


r=r (0) 

I  de 


0=0  r=e 

The  contribution  of  the  lower  limit  r=e  vanishes  in  the  limit  e-0. 
The  weight  functions  considered  in  this  appendix  are  admissible 
if  they  are  applied  to  residuals  which  are  sufficiently  smooth. 
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A 


Let  P  denote  the  weighted  residual  which  one  obtains  by 
applying  the  weight  function  -ty  to  the  potential  equation  (Eq.  (5)). 

p(x,y)  =  -  If  U-x,  n-y)  (<J>f  -  ( V2 )4>f2 )  +  (4>n)]d£dn 

S‘  '  '  (A20) 

i 

i 

The  region fi is  the  same  as  before  (Fig.  2).  The  weight  function 
is  regular  at  the  shock  (provided,  of  course,  that  the  point 
(x,y)  does  not  lie  at  the  shock) .  The  potential  equation  is 
supposed  to  be  satisfied  in  the  weak  sense,  in  particular  at  the 
shock.  The  shock  region  is  therefore  included  in  the  region  0 
and  no  discussion  of  the  shock  conditions  is  needed.  Carrying  out 
the  familiar  integrations  by  parts  one  obtains 

p(x,y)  =  -/  (£-x,  n-y)  [(<(>,.  (S,n)- ( V2)  4>f2(C,n)  )dn-4>  (£,n)d£]  \ 

30  c  c  C  n  ;  i 

+  f  h>  U  #  n)  (£-x,  n-y)  dn-tK  ( C-x,  n-y)  dU  (A21) 

an  Cn  : 

-  (  V2 )//  (C-x,  n-y)  <J>£ -  J  U,  n)dCdn 

With  the  use  of  the  definition  Eq.  (A9)  ,  the  second  integral  is 
transformed  as  follows 

/  4>U,n)  lip^^  (C-x,n-y)dn-i|<^n(^-x,n-y)dC 

d£2  ; 

=  -  /  <t>U,n)  [*  (S-x,n-y)dn  +  ^rn (t-x, n-y) df  (A22) 

30  ] 

D  ,  '  1 

=  -  /  <M£,n)-£  ( ihn ) dn  =  - 4>  ( C # n )  I  +  /  (£-x,n-y)  (4>  dn  +  4>  dfj 

an  an  n  1  I  ?  '  1  n  t,  ,  j 

limits  !  ; 

i  : 

The  portions  of  the  contour  30  ,  30o  and  30  form  a  closed  curve.  f  ■ 

■L  3  I 

The  contributions  at  the  limits  therefore  vanish.  The  same  holds  ' 

for  the  small  circle  30^.  Across  30  3  (portion  of  the  1  , 

wake)  ,  ^r=u  and  4>  =v  are  continuous ,  the  integral  therefore  vanishes 

t  n  1 

l 
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for  3Sl3.  The  first  two  integrals  in  Eq .  (A21)  evaluated  for  3Q2 

therefore  give 

/  -(u(£,n)  -  (  V2)  u2  (£,  n)  ]  ipf  (£-x,  n-y)  dn  +  u(£,n)ij>  (£-x,n  -y)d£ 

3  SI 

2 

+  /  v(£,n)[<Jj  (£-x,n-y)d£  +  ^  (r-x,n-y)dn]  (A23) 

Dfi2 

The  contour  3S14  (small  circle  around  the  point  (£,13)  *=  (x,y)) 
yields 

-  Ji^  (£-x,n-y)  [  (<j>^  (£,n-  (  Vz)  4^2  (£,n) )  <dn- 4>n  (£,n)  dn] 


r  ^cosG [u ( £ , n) - (  V2 ) u2 ( £ , n) rcos0  +  v ( £ , n) rsinO ] d0 


=  it  [  u  ( x ,  y )  -  (  Va)  u2  (x,y)  ]  (A24) 

Evaluating  Eq.  (A22)  for  one  obtains 
J>n<£-x,n-y)  [<j>ndn  +  <^d£] 

=  liro  /  r  ^sin0[  v( £, n) rcos0-u (£, n) rsinO ]d0 
£->■0 

0  =  0 
r=e 

=  tt  u(x,y)  (A25) 

Substituting  Eqs.  (A23) ,  (A24)  and  (A25)  into  Eq.  (A21)  and  com¬ 
paring  with  Eq.  (A17)  one  finds  that,  indeed, 

P(x,y)  =  Qx(x,y)  (A26) 


2  TT 

=  -lim  / 

£->0  Q 

0=0 

r=£ 
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APPENDIX  V 


DERIVATION  OF  THE  EXPRESSION  Qx  BY  APPLICATION  OF  THE 
WEIGHT  FUNCTION  -  ip ^  TO  THE  FORMULATION  IN  TERMS  OF  U  AND  V 

The  derivation  to  be  shown  here  has  of  course  great  similarity 
to  the  one  of  Appendix  IV.  There  are  some  differences  in  the  de¬ 
tails.  One  additional  question  arises.  In  the  formulation  of  the 
problem  in  terms  of  the  potential  <j>  this  quantity  is  always  assumed 
to  be  continuous  through  the  shock.  To  ensure  this  continuity  in 
a  computation  in  terms  of  u  and  v  one  must  add  the  postulate  Eq. 

(15),  i.e. 


[u]  +  (dxs/dy)  +  [v]+  =  0  (A27) 

A  second  shock  condition  is  given  by  conservation  of  mass,  Eq. 

(A14) 

[u  -  (l/2)uV  -  [v]  +  (dxs/dy)  =  0  (A28) 

First  we  show  that  the  partial  differential  equations  Eqs .  (12) 
and  (13)  viz 

^(u  -  (l/2)u2)  +  ^v  =  0  (A29) 

u  —  v  =0  ( A30 ) 

y  x 

and  the  shock  conditions  (A27)  and  (A28)  give  the  expression  Qx 
if  one  applies  the  weight  function  -  to  the  equation  of  con¬ 
servation  of  mass  (A29) .  Subsequently  we  shall  ask  whether  the 
postulate  Qx(x,y)  =  0  (for  all  values  of  (x,y) )  ensures  that  Eqs. 

(A27)  through  (A29)  are  satisfied.  It  is  assumed  that  the  u  field 
and  the  v  field  are  always  connected  by  Eq.  (A30) . 

To  derive  the  first  result  we  apply  an  integration  by  parts 
to  the  following  expression,  in  which  the  integrand  is  formed  with 
Eq.  ( A29 )  \i 
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p(x,y)=  -  //  ^ (C-x,n-y) [gj(u(£,n) - (1/2) u2 (£ ,n) )  +  ~v]dfdn 


=  -  /  (C-xfn-y)  t(u(c,n)-(l/2)u2(C,n))dn  -  v(c,n)drj 

3ft  s 


+  //  [i^  (£-x,n-y)  (u(Cfn)-(l/2)u2(rJ#n))  +  ^^n(C-x,n-y)v(C/n)  ]d£dn 


One  has 


J I  [r.r(£-x,n-y)u(£,n) 


=  // t-'l'  (S-x#n-y)u(f;,n) 

ft  nn 


=  /  ^  (£-x,n-y)u(£,n)d£ 
3ft  n 


+  ip ( £-x) v(c, n) ] d^dn 
+  i^n(£-x,n-y)  v(£,n)  ]d£dn 


+  i|>n  (£>y,n-y)v(s,n)dn 


+  J  J  <J>  (£-x,n-y)  (u  ( C# n)  -  v  (£,n)d£dn 
ft  n  n  4 

The  double  integral  vanishes  because  of  Eq.  (A30) .  One  thus  obtains 

p  (x,y)  =  /  -(u(c,n)-(l/2)u2(C,n)  (C-x,n-y)dn  +  u(£,n)i|'  (£-x,n-y)d£ 

3ft  4 

+  j  v(£,n)  (£-x,n-y)d£  +  i|>  U-x,n-y)dnl 

3ft  4  1 

-  (1/2  fj  \p  f  (C-x,n-y)  u2  (£,  n)  d£dn  (A31) 

ft  ^ 

The  contribution  of  the  contour  3ft^  (large  circle  which  moves  to 
infinity)  vanishes, so  does  the  contribution  of  the  integral  along 
the  wake  (contour  3ft^)  because  of  the  continuity  of  u  and  v.  The 
contour  3ft^  (small  circle  around  the  point  (£,n)  =  (x,y))  gives 
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2  I!  ( u  ( x ,  y ) 


( 1/4 ) u2 ( x , y ) . 


With  these  results  one  shows  that  indeed 


The  shock  need  not  be  considered  as  part  of  the  boundary  of 
ft  because  we  use  the  concept  of  weak  equality  in  Eq .  (A29). 

Nevertheless  let  us  introduce  a  contour  3ft,.  which  is  formed 
by  the  two  borders  of  a  cut  extending  along  the  shock.  (Traveling 
along  the  contour  one  must  leave  the  interior  of  ft  to  the  left.) 

One  therefore  travel,  in  the  positive  n  direction  on  the  upstream 
side  of  the  shock.  This  contour  gives  a  contribution 

/  'Pr  U-x,n-y)  [-(u(C,n)-{l/2)u2(fwn) )  +  v(£,n)  <d£/dn)  ]  +_  dn 

3ft  +  ^  s 

5 

+  /  Ij;  (£-x,n-y)  [u(E,n)  (d£/dn)  +  v(C,n)]+dn  (A32) 

3  ftt  n  S 

where  3ft^  is  the  upstream  portion  of  the  contour  9ft^.  The  two 
integrals  vanish  because  of  Eqs.  (A27)  and  (A28) . 

Next  we  assume  that  we  have  a  u  field  and  a  v  field  connected 
to  it  by  Eq .  (A30),  and  that  the  expression  P(x,y)  vanishes  for  all 

values  of  (x,y).  Each  value  of  P  for  a  given  (x,y)  arises  by  the 
use  of  a  weight  function  (£>x,  g-y)  .  By  an  argument  similar  to 
that  shown  at  the  beginning  of  Section  III  one  can  form  linear  com¬ 
binations  of  such  weight  functions  which  represent  weight  functions 
w(x,y)  with  finite  support.  We  choose  for  this  purpose  weight  func¬ 
tions  which  do  not  straddle  the  shock.  Then  it  follows  that  Eq . 

(A29)  is  satisfied  everywhere  (in  the  weak  sense)  in  the  region  where 
the  shock  is  excluded.  Therefore  the  residual  expression  (A31)  for 
the  region  ft  with  the  shock  excluded  vanishes  identically.  Let  us 
denote  this  expression  by  P(x,y).  P(x,y)  is  identical  with  P(x,y) 
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except  that  the  contribution  of  the  shock  boundary  must  be  taken 
into  account 


P(x,y)  =  P  ( x ,  y )  +  f  'I',  (P-x,n-y)  [-(u(f  ,n)-(l/2)u2(f.,n)  ) 

30 
‘  5 

+  v(c,,n)  (dts/dnl^  dn 

+  /  ,  h1  (C-x,n-y)  [ u ( f, ,  n )  (d; ;  /dn)  +  v(f,n)]+dn 

9fi5  1  s 


We  have  postulated  that  P(x,y)  =  0  and  we  have  shown  that  P(x,y)  =  0.  / 

The  sum  of  the  two  integrals  therefore  vanishes  identically  in 
(x,y) .  The  sum  of  the  two  integrals  represents  solutions  of  the 
Laplace  equation,  in  the  whole  x,y  plane  except  at  f°r 

-  ( :,-x,  n-y)  and  (  Z-x ,  n-y)  are  solutions  of  the  Laplace  equation. 

At  the  shock  this  function  has  jumps  in  the  normal  and  in  the 
tangential  derivatives  which  can  be  expressed  in  terms  of  the 
factors  of  ^ ^  and  .  But  the  above  mentioned  analytic  function 
is  identically  equal  to  zero.  The  jumps  are  therefore  zero  and 
these  two  factors  will  vanish.  Hence  it  follows  that  Eqs.  (A27) 
and  (A28)  are  satisfied. 
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APPENDIX  VI 

DERIVATION  OF  THE  EXPRESSION  Q FROM  THE 
PARTIAL  DIFFERENTIAL  EQUATION  FOR  u 

In  order  to  recognize  the  effect  of  the  nonlinear  term  we 
introduce 

q  =  (  Vz)  (u2)  (A3 3) 

This  quantity  can  be  regarded  as  a  source  term  in  the  Laplace  equa¬ 
tion.  The  formulation  of  the  problem  in  terms  of  u  is  given  by 
Egs .  (15),  (16),  and  (17).  Then  one  has 

u  +  u  -  q  =  0  (A34) 

xx  yy  x 

[u]+  -  [u2]*(cos23)/2  =0  (A3 5) 

and  as  a  substitute  for  Eq .  (17)  which  is  a  consequence  of  Eq .  (14) 

- [u]  +  tg$  +  [  v]  *  =  0  (A36) 

We  had,  Eq.  (7) 

tg3  =  -dxg/  dy 

In  Eq.  (A36)  it  is  assumed  that  one  has  constructed  a  v  field, 
which  is  connected  with  the  u  field  by  the  irrotationality  condition 

u  -  v  =  0  (A37) 

y  x 

in  the  manner  described  at  the  end  of  Section  II. 

The  difference  between  the  present  formulation  and  that  of 
Appendix  V  lies  in  the  fact  that  Eq.  (A37)  is  used  merely  to  deter¬ 
mine  the  v  field,  while  the  basic  problem  is  solved  in  terms  of  u. 

If  there  is  a  residual  in  Eq.  (A34) ,  then  there  will  also  be  a 


residual  in  the  expression  from  which  this  equation  is  derived. 
Let  the  latter  residual  be  r(x,y) .  Accordingly 


Ux  +  Vy  "  q  =  r(x'y) 


(A38) 


Then  one  obtains  for  the  residual  in  Eq.  (A34) 


u  +  u  -  q  =  r 

XX  yy  X 


(A39) 


The  expression  (A34)  is  not  the  equation  of  conservation  of  mass, 
but  only  its  x  derivative,  in  spite  of  its  form  it  is  not  a  conser¬ 
vation  law,  and  therefore  need  not  to  be  satisfied  as  one  traverses 
a  shock.  Test  functions  of  the  kind  used  in  the  theory  of  gener¬ 
alized  functions  cannot  be  used  if  their  support  region  straddles 
the  shock.  The  functions  tp  allows  one  to  construct  inadmissible 
as  well  as  admissible  test  functions.  We  denote  the  residual 
expressions  which  arises  from  Eq.  (A39)  by  P. 


P(x,y)  =  IS  iMS-x, n-y)r  (£,n)d£dn 

ft  ^ 

=  II  i^(C-x,n-y)  (u  (£,n)  +  u  (£,n)  -  q r(£,n))d£dn 

ft  11  n  ^ 


=  /  iHt-x,n-y)  (u  (f,,n)dn  -  u  U,n)df,} 

Oft  t,  n 


-  /  u(£,n)  (C-x,n-y)dn  -  (£-x,n-y)dU 

3ft  ^  1 


-  //  qf  (£, n)  ip  (£-x, n-y)  d£dn 

ft  ^ 


(A40) 


Here  the  contour  3ft  includes  the  two  borders  of  the  cut  along  the 
shock  (denoted  by  3fts) .  In  the  first  integral  on  the  right  we 
replace  u  by  v  by  means  of  Eqs.  (A37)  and  (A38) 
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/  iM£-x,n-y)  (u  (f.,n)dn-u  ('.,n)}d£ 

9ft  1  n 

=  -  /  iM£-x,n-y) (v  (£,n)dn  +  vr(C-n)dC  +  /  ip (£-x,n-y)q (C,n)dn 

1 1  S 

311  9ft 

+  j  i^(C-x,n-y)r(C»n)dn 

9ft 

=  /  vg,  n)  (ip.  (C-x,  n-y)d£  +  ip  (£-x,  n-y) dn)  +  /  <P  (£-x,n-y)q  (5, n)  dn 

3ft  9ft 

+  /  tp(£-x,n-y)  r  (£,n)dn 

3ft 

In  the  last  step  an  integration  by  parts  along  the  contour  9ft  has 
been  carried  out.  The  term  outside  of  the  integral  vanishes 
because  the  contour  is  closed.  The  integrals  in  Eq.  (A40)  over 
the  contour  at  infinity  9ft^  vanish.  So  does  the  integral  over 
the  wake  9ft 3  because  of  the  continuity  of  u  and  v  through  the  wake 
in  a  steady  plane  flow.  The  integral  over  the  small  circle  around 
the  point  (x,y)  is  treated  in  a  familiar  manner  and  one  obtains 
2 tt  u(x,y)  .  One  thus  obtains  from  Eq.  (A40) 

P  (x,y)  =  2Ttu(x,y)  +  /  v  (  £ ,  n)  ( *Pf  ( £-x,  n-y )  d£  +  <p  (  £-x ,  n-y )  d  n}d£  } 

3  2  3  ^  T> 

+  /  u(£,n) {-ip^  (£-x,n-y)dn  +  ip^ (£-x,n-y)d£) 

3  £2  2  3  ^  5 

+  /  iM£-x,n-y)q(C/n)dn  (A41) 

9ft  2  +  9ft  5 

+  /  ip  (£-x,n-y)  r(£,n)dn 


-  ff  qr  (£,n)il>  (£-x,n-y)d£dn 

ft  ^ 
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The  function  q  which  is  the  expression  for  the  nonlinear  terms 

enters  the  last  equation,  once  because  of  the  source  terms  in 

Eq.  (A34) ,  a  second  time  when  the  derivatives  of  u  are  replaced 

by  the  derivatives  of v  (by  means  of  Eqs.  (A37)  and  (A38)).  In 

the  original  differential  equation  the  nonlinear  terms  appear  in 

2  2  2 

the  form  q^  =  -(1/2)3  (u  )/3x'  ;  they  are  independent  of  the  local 
value  of  u  or  its  first  derivative.  The  rewriting  of  the  u 
derivatives  generates  a  second  nonlinear  term.  The  derivation 
has  been  carried  out  in  the  form  just  shown  to  make  clear  that 
the  term  so  obtained  is  not  solely  the  source  term  in  Eq.  (A34) . 
These  terms  can  be  combined 

-ff  q  (f,,n) 4* (?>x,n-y)d£dn  +  /  q(S»n)<M£-x,n-y)dC 

S3  f  3S3 

=  ff  q(£ ,  n)  !jv  (f-x,  n-y)  d£dn 

S3 

=  (1/2)  ff  (u2(£.,rH f(5-x,n-y)dedn 

S3  ”  ’  ^ 

This  expression  can  be  transformed  by  a  further  integration  by 
parts,  so  that  derivatives  of  the  unknown  function  u  no  longer 
appear.  In  evaluating  the  result  one  must  remember  that  one  deals 
with  the  region  S3  punctured  at  the  point  (x,y)  .  One  obtains 

+  (l/2)  ff  ly  u2  (f,,n)  (S-x,n-y)d£dr) 

S3  ’ 

=  (l/2)/u2  (f  ,n)  iK  U-x,n-y)do  -  (1/2)  ffu2  ( f  ,n)  4>r r  ( f-x, n-y)  d^dr, 

^  S3 

Here  the  double  integral  must  be  evaluated  for  the  punctured 

region  S3  (if  the  point  (x,y)  lies  within  S3).  The  last  form  has 

the  advantage  that  it  automatically  takes  6  functions  into 

account  which  might  for  instance  arise  on  the  left  side  of  the 

2 

last  expression  if  in  some  discretization  procedure  u  is 
approximated  regionwise  by  constants. 

The  discussion  just  carried  out  has  the  purpose  to  explain 
that  in  the  final  form  of  the  expression  the  nonlinear  terms  enter 
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in  a  rather  complicated  manner.  The  final  expression  for  P  is 
best  obtained  by  substituting  immediately  the  form  of  q. 


i 


P  (x,y)  =  ff  iMC-x,n-y)  {— ^5—  (u  U ,  n)  -  (1/2)  u2  (£ ,  n) )  +  u  (£,n))d£dn 

r\  ^  r  **  '  I U 


=  f  iMC-x,n-y){  tI-  (u(C,n)-(l/2)u2(C,n))dr1  -  un  (C,n)dO 
3ft  n 

(A42) 

-  /  { (u(£,n)-(i/2)u2(£,n) )  (<K  (C-x,y-y)dn  -u (£ , n) (C-x, n-y) d£ } 

3ft  t,  n 


~ (1/2)  ff  *  (€-x,n-y)u2(£,n)d£dn 

ft  ^ 

In  the  first  integral  on  the  right,  u  is  replaced  by  v  by  means 
of  Eqs.  (A37)  and  (A38) .  (A  corresponding  step  has  been  carried 
out  before.)  One  obtains 


-  /  iMS-x,n-y){v  U,n)dn  +  v_(£,n)dU+  /  (MS-x,n-y)  r  (£,n)d 
3ft  n  ^  3ft 


=  f  v(c,n)  {<Pr  (^-x,n-y)dC  +  <Pn  (£-x,n-y)dn>+  /  <M€-x,n-y)r(£;,n)dn 
3ft  ^  n  3ft 


Again  tne  terms  outside  of  the  integral  vanish  because  the  contour 
3ft  is  closed.  The  integral  for  the  small  circle  around  the 
singular  point  (x,y)  is  treated  as  above.  One  obtains 

P(x,y)  =  a  [tr  (u(x,y)  -  (1/4)  u2  (x,y))  ] 


+ 


/ 

3ft2  +  3ft^ 
/ 

3ft2+3ft^ 


^ (C-x,  n-y)  [  (u(£,  n)  -  (l/2)u2  (£,  n)  )dn-v(£,  n)d£)  (A43) 


if,n  (£-x, n-y)  fu(f;,n)dC+v(^,n)dn)  ] 


(Continued 
next  page.) 
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-  (1/2)  ff  ip  (C-x, n-y)  u  (£ , n)d£dn  +  /  iMf-x,n-y)r  (f;,n)dn  (A43) 


€€ 


3Q 


(Concluded) 


with  a=  2  for  (x,y)  within  ft 

a=  0  for  (x,y)  outside  ft. 

The  contributions  of  the  shock  contour  3ft_  vanish  because  of  the 

D 

shock  conditions  (A35)  and  (A36) .  Therefore,  by  comparing  with 
Eq.  (A17)  and  using  Eq.  (A26) 


P(x,y)  =  Q  (x,y)  =P(x,y)  -  /  iJj  ( £-x,  n-y )  r  ( g ,  n)  dn 

x  3ft 


This  result  is  to  be  expected  on  the  basis  of  the  following 
reasoning.  Using  the  original  definitions  one  has,  from  Eqs.  (A20) 
and  (A38) 

P  (x,y)  =  ff  ip  (Z-x,r\-y)rf.  (f^nJdf.dn 

ft  ^ 

=  /  ip  (C,n)r(C,n)dri  -  ff  r(£,n)iMS-x,n-y)d£dri 

dfi  ft  ^ 

=  /  iMS-x,n-y)r(£;,n)dn  +  P(x,y) 

3ft 

The  residual  expressions  P  (arising  from  a  formulation  in  terms  of 
u)  and  P  (arising  from  a  formulation  in  terms  of  <j>  or  in  terms 
of  u  and  v)  are  closely  related  to  each  other  but  not  identical. 

This  is  caused  by  the  fact  that  Eq.  (A39)  with  rx=0  merely  guarantees 
that  the  function  r  which  appears  in  Eq.  (A38)  is  a  function  of  y 
on^y*  !n  working  with  Eq.  (A39)  we  impose  the  additional  require¬ 
ment  r  =  0  at  infinity  and  use  the  original  expression  I. 

To  treat  the  limiting  case,  in  which  the  point  (x,y)  approaches 
the  contour  3ft2  from  the  interior  of  ft  we  rewrite  the  integrals 
over  9^2  in  Eq.  (A43)  as  follows. 


(A44) 


/  u(C,n)  {-4>r  (Z-x,  n-y)dn  +<l>  (£-x,n-y)d£} 

3«2  5  " 


+  f  v(C , n)  (C-x, n-y) d£  +4*  (£-x, n-y)  dn) 
dQ2 


+  (1/2)  /  u2(C/n)^F(?-x,n-y)dn 


Let  3  be  the  angle  of  the  *nner  normal  of  the  region  ft  with 
the  x  axis  (Fig.  5) .  Then,  for  the  unit  vector,  in  the  direction 
of  the  inner  normal 

e  =  e cos  8  +  e„sin3 
n  x  y 

and  for  the  unit  vector  in  the  direction  of  the  tangent  to 
pointing  in  the  direction  in  which  the  integration  proceeds 

eg  =  exsin6  -  e^,cos3. 

If  do  is  the  line  element  of  8ft  one  has 

~4§  =  sin  6,  ^  =  -cos0 


then 

a5r  cos3  +  VinS  =  air  +%if 
3tT  =  VinB  -*n  cose  =  +^n  la 

^  =  fnT  CosB  +  3t  SinB 

i?  =  TrT  Cos2b  -  it"  SinBcose* 


(A45) 


} 

I 


) 


(A46)  I 

\ 

i 


The  expression  (A44)  then  assumes  the  form 


/  [u(S,n)  -  (l/2)cos20  (o)  u2  (S  ,n)Mdi^/dn)do 

3^2 

+  /  [v{£,  n)  -  (1/2)  sing  (a)  cos3(o)  u2  (C  ,  n)]  (di^/dt)  do 

3n2 


(A47) 


We  defined 

i|*  =  log  r(£-x,  n-y)  . 


We  consider  the  limiting  process  where  the  point  (x,y)  approaches 
a  point  (XC'YC)  along  the  inner  normal.  Applying  a  procedure 
familiar  from  potential  theory,  one  treats  the  integral  along 
3£2_  in  a  small  vicinity  of  the  point  (x  ,y  )  separately.  In  the 
limit  where  (x,y)  coincides  with  Xc/Yc)  the  integrand  of  the  first 
integral  in  Eq.  (A47)  remains  bounded  in  the  outer  region  (where 
the  point  (xc,yc)  is  excluded.  For  the  small  region  which 
contains  the  point  (XC»YC)  one  obtains 

/(diJi/dn)do  =  -it. 

The  second  integral  in  Eq.  t,A47)  must  be  evaluated  as  Cauchy 
principal  value. 

One  thus  obtains,  from  Eq.  (A17) ,  for  the  point  (x(s),y(s)) 
on  the  contour 

P(x,y)  =rr(u(x,y)  -  (l/2)sin23(s)u2(x,y)  (A48) 

+  {u(c,n)-(l/2)cos2B(s)u2(£,n)  H-i|v  (c-x,n-y)dn+i|»,  (?-x,n-y)dO 

o«2  ^  ^ 

+  /  ( v (r , n)  -  (1/2) sing (s) cosB  (s)  u2  ( C , n)  (S-x,n-y)d?+i|>  (5-x,n-y)df,) 

-  (1/2)  If  (^-x,n-y)  u2  (S,n)dCdn 

n  ’ 
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To  avoid  an  ambiguity  in  the  signs ,t|>^  and  i^,  instead  of  and  1 J»t 
have  been  reintroduced,  although  the  and  express  the  nature 
of  the  integrands  more  clearly. 

The  evaluation  of  the  double  integral  requires  some 

precautions  because  of  the  singularity  introduced  by  the  factor 

in  the  integrand.  The  region  under  consideration  is  the 

punctured  region  ft.  In  a  specific  case  the  procedure  depends 

2 

upon  the  chosen  representation  for  u  .  To  show  what  will  happen 

2 

we  consider  the  case  where  u  is  constant  within  an  element 
bounded  by  straight  lines.  Then  one  obtains 


- ( 1/2 ) u2 


//  i d£dn  =  -(l/2)u2 

n' 


/ 

a  q 


y^dn 


where  Q '  is  the  element  under  consideration  and  its  boundary. 

If  the  point  (x,y)  lies  within  the  element,  then  for  the  small 
circle  around  this  point 


/  y^dn  =  -it. 

For  the  remainder  of  the  contour  it  does  not  matter  whether  the 
point  (x,y)  lies  inside  or  outside  of  the  contour. 

We  introduce  n~y  =  n'.  Then  one  has 

=  - i! - 2 

Let  the  contour  between  the  points  A  and  B  in  Fig.  6  be  given  by 


=  a  -  n'tgS 


where  3  is  the  angle  between  the  inner  normal  to  the  boundary  AP 
and  the  x  axis  (Fig.  6.)  Then  one  has  to  evaluate 


[•IT 


(a-n'tgg)  +n 


n 

A 


5  !-n  n* 


The  direction  of  integration  is  chosen  in  such  a  manner  that 
the  interior  is  to  the  left  as  one  travels  along  the  contour.  Thi 
is  an  elementary  integral. 


mo  (— — ^  -a  sin3)  (^-7) 

/ - — HSBB - dq'  =  -singcosB  /  — -? - 2 — 2  2~ 

rt?  (a-n’tg  )2+n'?  nl  (5-^3  -asinB)  +a  c'os  3 


+cos^  3  / 


a  cos  6 


nA  t(^i6  -asin3)/acos3]  +1 


=  -  sin  3cos3  (1/2)  logf  ”  asing)  2+a2  cos23)  j 


n  .  n  1  b 

•>  — — a  -a  sm3 

+  cos23  arctg  - 5 - 

J  a  cos6 


One  has,  either  from  Fig.  7  or  by  direct  computation, 
r 2  =  n'2  +  K'2=  ( — \  -  a  sin3) 2  +  a2  23 


0A=  arctg 


j-  a  sin 3 
a  cosB 


Thus,  finally  as  the  contribution  of  the  boundary  AB 


- (1/2) u2 [-sin3  cosB  log  —  +cos23(0  -  0  )] 

rA  B  A 


We  apply  this  result  to  several  special  cases. 


Let  the  point  (x,y)  lie  in  the  middle  of  a  square  (Fig.  8a). 
One  then  obtains 

-  ( 1/2 )  u2  //  i|»  dCdn  =  -  \  u2[-tt  +tt]  =  0 

a  ^ 

The  first  factor  comes  from  the  contour  around  the  point  (x,y) . 

The  second  one  from  the  contour  AB  and  CD.  Incidentally,  one 
sees  that  the  exclusion  of  the  point  (x,y)  by  a  square  is 
equivalent  to  the  exclusion  by  a  small  circle. 

For  a  narrow  rectangle  (extending  in  the  y  direction, 

Fig.  8b)  one  obtains 

(-1/2)  u2  //  dCdn  =  -  \  (-7T  +  2tt)  =~j~  ■ 

If  the  point  (x,y)  approaches  the  contour  AB,  Fig.  8c,  then 
the  contribution  of  the  contour  AB  cancels  the  contribution  of 
the  small  circle.  In  this  manner  one  can  discuss  the  problem 
and  contours  of  different  kind;  for  instance,  the  contribution 
of  an  oblique  contour. 

One  finds  that  the  local  contribution  of  the  double  integral 
in  Eq.  (A47)  will  cancel  the  local  contribution  (-l/2)sin28  u(x,y) 
and  also  the  local  contributions  of  the  second  contour  integral 
(in  which  the  Cauchy  principal  value  is  to  be  formed) . 


APPENDIX  VII 


RELATION  BETWEEN  THE  SHOCK  CONDITIONS  EQS .  (14) 

AND  (15)  AND  LOCAL  EXPRESSIONS  Q  or  Q 

x  y 

It  has  already  been  shown  in  Appendix  V  that  the  shock  con¬ 
ditions  (14)  and  (16)  are  satisfied  if  Qx(x,y)  =  0  identically  in 
x  and  y.  In  this  appendix  it  will  be  shown  that  if  Qx  =  0  or 
0^  =  0  locally  at  points  immediately  upstream  and  downstream  of 
the  shock  then,  respectively,  Eq.  (14)  and  a  linear  combination  of 
Eqs.  (14)  and  (16)  will  be  satisfied.  In  other  words,  if  one  uses 

in  the  computation  either  conditions  for  Q  or  Q  then  one 

x  y 

relation  derived  from  the  two  shock  conditions  is  satisfied  even 

if  Q  =  0  or  Q  =0  only  at  points  immediately  upstream  and  down- 
x  y 

stream  of  the  shock. 

This  discussion  is  primarily  of  interest  for  the  classical 
use  of  the  integral  equation  method.  Accordingly,  it  suffices  if 
we  restrict  it  to  the  specialized  equations  (25)  and  (26)  although 
the  results  have  general  validity.  Let  us  first  consider  the 
condition  derived  from  Q^.  We  evaluate  the  expression  (26)  for 
points  (x,y)  approaching  the  shock  from  the  upstream  and  the  down¬ 
stream  side  along  a  line  y  =  const.  Let  the  point  of  the  shock 

reached  in  the  limit  be  x  (y) .  A  discontinuity  in  the  expression 

s  2 

Q  can  be  caused  by  the  terms  2uu(x,y),  —  ( tt/2 )  u  (x,y)  and 

X 

-(1/2)  ff  u2(C,  n)tr  (C-x,n-y)d^dn 

As  always  the  region  Q  is  part  of  the  £,n  plane  punctured  at 
U,n)  =  (x,y).  We  enclose  the  point  (xg,y)  by  a  small  box  (Fig.  9). 
In  the  limiting  processes  where  one  approaches  this  point  from  the 
upstream  and  downstream  side  the  singular  point  (x,y)  moves  with¬ 
in  the  box.  We  consider  the  limits  of  the  double  integral  for 

x  =  lim  (x  -e)  and  x  =  lim  (x+e);  e>0 
t-0  S  e->0 
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a 


In  this  limiting  process  the  contributions  of  the  double  integral 
from  outside  of  the  box  will  be  continuous  as  x  moves  through  the 
shock.  Only  the  region  within  the  box  need  to  be  considered.  Let 
ft  +  and  ft  be  the  portions  of  ft  within  the  box  upstream  and  down¬ 
stream  of  the  shock.  Let  u+  and  u  be  the  values  of  u  upstream 
and  downstream  of  the  shock.  The  difference  of  the  double  integral 
for  the  area  ft+  +  ft  with  x  once  approaching  the  shock  from  the 
upstream  side  and  a  second  time  from  the  downstream  side  is  given 
by 


D 


/ /  (u+)  n)  (  £-xs 

ft  + 


(y) +e , n-y) dCdn 


+  //( u  )*( £,n)  ij>  .  (£-xs  (y)  +e ,n-y)d£;dn 


,  ,  -  2 
-  //  (u  ) 

ft” 


(  C-xg (y) -e , n-y) d£dn 


In  evaluating  the  first  and  fourth  of  these  integrals  one  must 

take  the  presence  of  the  singular  points  at  x  -e,y  and  x  +e,y 

s  s  + 

into  account.  One  can  take  the  box  sufficiently  small  so  that  u 
and  u  can  be  considered  as  constant.  (one  can  indeed,  show  that 
higher  order  terms  in  the  development  of  u  do  not  contribute  to 
a  discontinuity.  Then  the  integration  with  respect  to  £  can  be 
carried  out  in  the  above  integrals  and  one  obtains  only  integrals 
over  the  contour.  The  evaluation  of  the  contours  around  the 
singular  points  is  familiar.  One  obtains 

-tt  ( u  )  and  -it  (u  )  . 

Therefore 
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12 


■t 


D  =  -ir(u  +  )  2  +  (u  +  )2  /  4>r(^-x  (y)  +r  n-y)dn 


3ft 


+  ' 


+  (u  )2  j  ipf(^-xs(y)  +c,n-y)dn 


3ft 


_  2  _  2 

-  (u  )  /  ijj  j.  ( £-xs  (y )  -e,n-y)dn  +  tt(u  ) 

dn~ 

Here  3ft  +  and  3,0  denote  the  contours  of  the  areas  ft+  and  0  ,  but 
with  the  exclusion  of  the  contours  around  the  singular  points. 

These  contours  include  the  borders  of  the  cut  made  along  the  shock. 
At  the  outer  boundaries  of  the  regions  ft+  and  0  namely  £  =  x^, 

Yl  <  n  <  y2  and  c  =  x2;  Yl  <  n  <  y2  (pig.  9)  one  has  obviously 

lim  ip  (F,-x  (y)  +  e,n-y)  =  lim  ( £-x^  (y) -e  ,  n-y) 
e+0  **  s  e->-0  ^  s 


In  the  limit  e->-0  these  parts  of  the  contour  do  not  contribute  to 
the  difference  D;  only  the  integrals  along  the  shock  are  left. 

One  observes  that  this  integration  is  carried  out  in  opposite  direc¬ 
tions  for  the  contours  3ft  +  and  3ft  ,  on  the  upstream  side  in  the 
direction  of  increasing  n .  One  therefore  obtains 


D  = 


ii  t  \  2 

[ (u  )  - i 


u  )  2  ]  {  -it 


U-xs(y)+e,n-y)dn 


^2 

-  )  <J/£(S-x  (y)-e,n-y)dn> 

Yl  ’ 
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One  has 


i ( C-x, n-y) 


(  £~X) _ 

( C-x) 2+(n-y) 2 


We  set  for  the  evaluation  of  the  integrals 
£-xg  (y)  =  f, ' ,  n-y  =  n’ 

and  for  the  limits  of  the  integral 

nj  =  y-^-y,  n2  =  y2-y  • 

Notice  that 


ni  <  o  <  n2 

For  the  present  purpose  it  suffices  if  one  replaces  the  shock  by 
a  straight  line 

=  -n'  tg  3 


where  (as  in  Section  2)  3  is  the  local  angle  between  the  downstream 
pointing  normal  to  the  shock  and  the  x  axis.  Then  the  first 
of  the  above  integral  assumes  the  form 


I 


1 


_ +  £ _ 

(C'+e)2+n'2 


dr)  '=  / 

n-! 


-tg3n  +e 
(-tg3n'+e)  +n'2 


We  show  the  essential  steps  in  the  evaluation  of  this  elementary 
integral 
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(-tgBn  +e)2  +  n' 2  =  e.2cos23[( — ~  ■  2  ~  ~  tgB)2  +  1) 

ecoszB 


Setting 


— - 2 - tg(3  =  v 

GCOS  3 


one  obtains 


V2  2 
T  _  (  -sinBcosj3v+  cos  3  , 

“  J - 9 -  dv 

vx  1+V 


where  the  limits  of  the  integral  are  given  by 


v,  = 


ecos  3 


2 - tgB  and  v2  =  - ^ - tgB 


cc os  B 


For  e->-0,  e>0  one  has  v2-  +  03 

The  contribution  of  the  first  summand  in  the  integrand  of  1^ 
vanishes  because  of  antisymmetry  and  one  obtains 

T  2_ 

1^  =  cos  Bti 

The  second  integral  in  the  expression  D  is  treated  in  the  same 
manner,  but  >  is  replaced  by  -e .  Therefore  v^  =  +°°  and  v2  =  -°° 

One  thus  finds 

D  ~  «u+)2  -  (u  )2)(-ir  +  2ttcos2B) 
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■-■a* 


Forming  the  difference  of  the  expressions  for  Q  ,  (Eq.  25)  up¬ 
stream  and  downstream  of  the  shock  in  the  limit  e->-0,  one  obtains 
with  the  notation  of  Section  2 


2tt[u]+  -  tt/2  [  u2  3  ^  -  ( 1/2 )  D 


-  2  7T  [  U  ]  _  -  [u  ]_  ( 7T  /  2  -  tt/2  +  TTCOS  B)=  0 

This  is  the  condition  Eq.  (16) . 

An  analogous  procedure  carried  out  in  Eq.  (26)  requires  the 
evaluation  of 


-d/2)  //  u2(^,n)^Sn  U-x,n-y)d?dn 


After  an  integration  with  respect  to  K  one  arrives  again  at 
integrals  around  the  contours  3p+  and  .  No  contribution  due 
to  the  singular  point  (x,y)  appears.  One  arrives  at 

n2 

-  i[u2]*  {  j  t|>  U-x  (y)+e,n-y)dn 
n 


-  J  ip  (C-xs  (y) -s,  n-y)dn  } 
nl 


Here 


(C-x,n-y)  = 


(£-x)  +  (n-y) 


In  the  integral  one  makes  the  same  substitution  as  before, 
finally  obtains  from  Eq.  (26) 

2  [v]+  +  [u2]+  sinB  •  cosB  =  0 
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The  conditions  derived  for,  therefore  lead  to  a  linear  combination 
of  the  shock  conditions  (14)  and  (16) . 

The  full  shock  conditions  are  obtained  only  if  simultaneously 
Qx  =  0  and  =  0  for  points  immediately  upstream  and  downstream 
of  the  shock.  The  application  of  only  Qx  =  0  or  =  0  yields 
the  full  shock  conditions  only  if  the  respective  expression  van¬ 
ishes  globally- 


APPENDIX  VIII 


DERIVATION  PERTAINING  TO  EQUATIONS  (35)  AND  (36) 

The  evaluation  of  the  normal  derivative  of  <f>  in  Eq.  (35) 
requires  some  caution  because  the  singularity  which  arises  in  the 
derivatives  of  \p  as  the  point  (x,y)  approaches  the  contour  3Q. 
Therefore,  we  show  some  details.  Let  0(s)  be  the  angle  of  the 
inner  normal  to  the  contour  of  the  distant  field  with  the  x  axis 
(Fig.  5) .  In  the  contour  integrals  the  variable  s  will  be  replaced 
by  an  umbral  variable  a.  We  assume  that  s  (and  a)  increase  in 
a  direction  of  travel  for  which  the  distant  field  is  to  the  left. 
The  contour  is  described  by  x  =  x(s)  and  y  =  (s) .  We  define  unit 
vectors  in  the  direction  of  the  normal  and  in  the  tangential 
directions . 


Then  one  has 


ipf  (£  (o)  -x,  n  (a) -v  )  =  ij^cos  3  (o)  +  i^tsinB(o) 


i^n(C  (o)-x,n(o)-y )  =  ^nsin  3(a)  -  iptcosB(a) 


Now  from  Eq.  (34) 


4>  (x,y)  =  (2it)  -1  [-/  f  (a)  ip  (£  (a)  -x,  n  (a)  -y)  da-ffq  (£ ,  n)  <Pr  (£-x,  n-y)  d£dn 
x  3fi  ^  4 

+fr(x,y)  sin0(x,y)]  (A50) 


<p  (x,y)  =  ( 2  tt  )  —  1  [-/  f  ( o )  ip  (C(a)-x,n(a-y)da-//q(C,n)^r.  (C-x,n-y)d^dn 

y  3^  n  n  n 

-Tr(x,y)  ^  cosO(x,y)]  (A51) 


with  r2  =  x2+y2  sinO (x,y) =y/r ,  cos0(x,y)  =  x/r. 

The  above  expressions  for  ip f  and  are  substituted  in  the 
integrals  over  We  form 

c()  (x,yposS(s)  +  (x,y)sin3  (sM2tt) -1{_/f  (a)  (o,x,y)cos  (3  (o) -3  (s)  )do 
x  Y  dQ  n 


-  /  f ( a ) ^  (a,x,y)sin  ( 3 (a) -3 (s) ) do 

-cos3(s)  ff  crU , n)  Pc  (£-x,  n-y) d£dn 

ft 

-sin3(s  )ff  q(C»n)<l'n(C”x,n~y)d^dn 

ft  n 

-  ( r/r) sin ( 3 (s) -0 (x,y) ) } 


and  '  f 
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<f>x  (x,y)  sinB  (s)  -  <£y  (x,y)  cosB  (s) 


-1. 


=  (2 it)  i  -  /  f(o)ip  (o,x,y)sin  (B(s)-B(o)do 

9n 


-  /  f  (o)ip  (o,x,y)  (~cos  (B(s) -B  (o))  da 

an 


-  sinB  (s)  //  q(C/n)i|ir  (C-x,n-y)dCdn 


n 


+  cosB(s)  //  q(£,o)if)  (£-x,o-y)dCdn 


n 


-  (T/r)  cos  (8  (s)  -0  (x,y))} 


Now  one  carries  out  in  the  last  formulae  the  limiting  process 
where  (x,y)  approaches  a  point  of  the  contour  (x(s),y(s)).  The 


results  are  denoted  by  <(>n(s)  and  <J>^(s) 


.-1, 


I)  (s)  =  (f(s)/2)  +  (2rr)  {-  /  f  (cr)  ifi  (a  ,  (x(s)  ,y  (s)  )  Cos (g (a) -B (s) da 


9fi 


•  /  f(a)i//(a,x(s)  , y (s) )  sin  (8  (o) -3  (s))da 

an 


( A52 ) 


-  cosB  (s)  //q(C»n)  i>r  ( T, — x  ( s )  ,n-y(s)  )d£dn 


n 


sinB  (s)  ff  q(C»nH  (S-x(s)  ,  r|-y(s))d£dn 

n  n 


-  (T/r)  sin ( B (s)  -  0 (x (s) ,y  (s)  )  ) 


-1, 


(J>  (s)  =  (2tt)  {  -  /  f  (o)  i|j  (o,x(s)  ,y  (s)  )  sin  (  B  (o) -B  (s)  )  da 

u  an  1 


+  /  f  (a)  i/j  (o,x  (s)  ,y  (s)  )  cos  (3  (a) -B  (s) )  da 

an  c 


(A53) 


-sinB(s)//  q(Cf  n)  i(ir  (?-x,n-y)dCdn 


n 


(Continued  on  next  page.) 
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r 


+  cos 6(o)//  q(f;,n)ip  U-x,n-y)d£dn 

ft  n 

-(r/r)cos  (s(s)-e(x) (o) ,y  (s) )  } 

The  present  representation  of  4>n(s)  and  ^(s),  which  uses  i| >  and 
i l>t  may  appear  less  direct  than  a  representation  in  $  and  ^  but 
it  shows  more  clearly  the  limiting  form  of  the  equations.  If 
the  point  (o)  /  n  (cr) )  approaches  the  point  (x(s),y(s))  (in  other 
words,  as  o  approaches  s)  .  The  derivative  ipfc  behaves  as  (a-s)  ^ ; 
integrals  containing  must  be  evaluated  in  the  sense  of  Cauchy 
principal  values  (unless  there  is  another  factor  in  the  integrand 
which  tends  to  zero  in  the  same  manner).  Furthermore, 

lim  tp  =  tp  (o,x(s)  ,y  (s)  )=  -(2d3(s)/ds)  1 
s  n 

(One  will  remember  that  s  increases  in  a  direction  in  which  one  has 
the  distant  field  to  the  left.  d3/ds  is  positive  if  the  distant 
field  is  convex  at  the  point  s.) 

If  one  evaluates  the  double  integrals  in  the  form  used  in 

Eq.  (35)  then  it  does  not  matter  whether  the  plane  ft  is  punctured 

or  not.  If,  however,  one  transforms  it  by  integrations  by  part 

2 

so  that  the  derivatives  of  u  in  q  vanishes,  then  one  must  remember 
that  one  deals  with  the  punctured  plane  in  the  limit  where  the 
singular  point  approaches  the  boundary  from  the  inside  of  the 
distant  field. 

ffq (C,  n)  tpr  (C-x,  n-y)  d£dn=  (1/2  )//Tf  (u2  (£ ,  n))  :Pf  (C-x,  n-y)  dCdn 

ft  ft  ^ 

=  U/2)  f  u2  (C,n)  ipr  (£-x,n-y)dn  -  (1/2)  ffu2  (£,  n)  (S-x,  n-y)d£dn 

3ft  ’  ^ 

=  -(tt/4)u2  (X,y)  +  (1/2)  /  u2  (€  (o)  n  (o)\pf  (£  (o)  -x,n  (o)  -y)  dn  (o) 

3ft2  5 

-  (1/2)  ff  u2( » n)  r  ( r~x,  n-y)  df,dn 
ft  ’  ’ 

lim  x-+x(s) 

y^yis) 

We  have  found  above 


84 


<l>£  (/  (j)  -x,  n  (O)  -y)  =  i^ncosB(o)  +  iptsin3(o) 


and  dn/do  =  -cosB(o).  Then 

(1/2)/  u2(U°)  ,n(o))4>r  (€(a)  -x,  n  (o)  -y)  dn  (o) 

a«2 

=  d/2)  /  u2(f.  (a)  ,  t|  (n)  )  {('/  (o,x,y)cos(Uo) 
du2  1 

2  2 

+  i|j  (o,x,y)sinP(a)  ]  (-cosB(a)  )  dn  =  -(n/4)cos  H  ( s )  u  (x(s),v(s)) 
-(1/2)  /  u2  (/(a)  ,n(t>) )  ]'(■  (o,x(s)  ,y(s)  )cos2h(o)+  .|i  (o,x(s)  ,y(s)sinfi(o)cosii(a)  ]do 

an. 

The  portion  of  the  integral  involving  ip  is  evaluated  as  Cauchy 
principal  value.  For  the  evaluation  of  the  double  integral,  see 
Appendix  VI. 

In  a  corresponding  manner  one  obtains 
//q(/,n)^n(f-x,n-y)df;dn  =  ^  //^  u2  (r. , n)  ipn (f,-x,n-y) d/dn 

=  i-  /  u2  (c,n)  (C-x,  n-y)dn  -d/2)//  u2  (f,,n)^fn(C-x,n-y)d/dn 

2  afi  n  J 

=  i-  /  u2  (/ ,  n)  (ip  sin(3  (a)  -ip  cos3  (a)  )  (-cos(3  (a)  )  do 

an2 

-  (1/2)//  u2(f.,n)'J'Cn  (f“X, f-y) d/dn 

=  +  ( n/4)u2(x  (s)  ,y  (s))sin  0(s)cos3(s) 


+  \  /  u2  (/ (a)  ,n  (a)  )  [-ip  (a,x  (s)  ,y  (s)  )  sing  (a)  cos6  (o) 

3«2 

+  <('t  (|> ,  X  (s  )  ,  V  (s)  eos2rt  (o  )  )  do 


-(1/2)  //  u2  (C,nW’  ,  (f-x(s)  , n-y  (s)  )ar,dn 

a  n 

limx+x  (s) 
y=y (s) 


f 

i 

I 

i 

In  practice  one  will  use  a  contour  which  consists  of 
straight  lines  oriented  in  the  x  ana  y  directions  and  considerable 
simplifications  will  be  encountered,  because  cos8  and  sin  3  will 
then  be  either  zero  or  one. 


J 


7k 


1 


I 
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Figure  3.  Notation  Used  in  Demonstrating  the 
Continuity  of  Q  as  the  Point  (x,y) 
Passes  Through  a  Shock. 
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UNCLASSIFIED 

2  0.2  I 


Figure  4 


Figure  5 


Coordinate  Systems  Used  for  the  Differentiation 
of  a  Double  Integral. 


Limiting  Case  Where  a  Point  (x,y)  of  the  Interior 
Approaches  a  Point  (x  ,yc)  of  the  Contour. 


Figure  6.  Notation  for  the  Evaluation  of  a  Double  Integral. 


Geometric  Interpretation  of 


Figure  11.  Derivation  of  Conditions  Along  9 ft-  if  Sources 
are  Present  in  the  Moderately  Distant  Field. 
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